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EXECUTIVE SUMMARY 


The Screem-J4D is a remotely piloted airplane that has been designed to fly at 
a chord Reynold's number of 100,000 while performing figure-8 maneuvers 
in a restricted area. It has a high aspect ratio main wing with a conventional 
empennage giving it a "sailplane" appearance. A specifications table and 
three-view drawing accompany this summary. 

The flight plan calls for ascent to cruise altitude at 20 ft and then perform 
three figure-8 turns around pylons. These pylons will be separated by fifty 
yards and stationed inside a sports facility. Once completed, the pilot is to 
make use of any remaining power by loitering before landing the plane. 

The propulsion system of the J4D consists of a propeller-electric motor 
combination with the prop mounted at the front of the fuselage. The 
Tornado 10 inch diameter, 6 inch pitch, two-bladed propeller is powered by an 
Astro 05 electric engine with 7 AA Nickel-Cadmium batteries. The system is 
capable of 50 watts of power and has throttling capabilities. Of the propellers 
with available data, the 10-6 was best fitted for the take-off distance and 
maximum current draw constraints. The 05 engine was chosen for being 
most lightweight while still supplying adequate power. 

In order to provide sufficient lift for low speed travel, the J4D has an aspect 
ratio of 11.72 with an 8.2 inch mean chord. The wing consists of a spar and rib 
construction with Micafilm skin. Its low-mount and dihedral in combination 
with the vertical tail were prescribed to allow for excellent maneuverability. 
A major problem, however, is that due to the nature of a low-speed mission, 
there is little margin for error between the cruise and stall velocities. As far as 
center of gravity travel is concerned, a square fuselage will internally contain 
the servos, engine, and so that the CG of the airplane is kept at about 33% of 
the chord. 


A combination of directional and longitudinal control will enable the J4D to 
perform the figure-8 maneuvers. However, in order to avoid the 
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construction and servo weight of ailerons, the rudder was designed to be over 
one-half the size of the vertical tail to insure that the proper roll motion 
could be attained. 

After a semester of hard work and discussion, the members of Group D feel 
that this proposed is a well-established answer to the station keeping mission. 
We are confident in the success of our design and believe that the Screem-J4D 
is worthy of continued support. 


Aerospace Engineers for the Twenty-First Century 
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SCREEM - J4D SPECIFICATIONS 


PERFORMANCE 


Endurance 

3.9 min @ 23 ft/s 

Range 

5500 ft @ 23 ft/s 

Stall Velocity 

18.9 ft/s 

Maximum Velocity 

58.6 ft/s 

Maximum Rate of Climb 

6.3 ft/s 

Maximum Power Available 

50 Watts 

Maximum Bank Angle 

30 degrees 

Maximum Roll Rate 

28.8 degrees /s 

Minimum Turning Radius 

28.6 ft 

Minimum Glide Angle 

-4.1 degrees 

Minimum Takeoff Distance 

42 ft 

Maximum Cl 

1.3 

Maximum Angle of Attack 

14.5 degrees 

CONFIGURATION 


Wing Span 

96 in (8 ft) 

RPV Length 

37 in (3.1 ft) 

Weight 

3 lbs (48 oz) 

Wing 

Airfoil Section 

NACA 4415 

Angle of Incidence 

3.5 degrees 

Aspect Ratio 

11.72 

Chord 

8.2 in 

Dihedral 

10 degrees 

Taper Ratio 

1.0 

Fuselage 

Cross section 

4 m X 4 in 

Payload Volume 

25.17 cu in 

Empennage Airfoil Sections 

Flat Plate 

Horizontal Tail Area 

.63 sq ft 

Elevator Area 

.21 sqft 

Vertical Tail Area 

.38 sq ft 

Rudder Area 

.21 sq ft 

Motor 

Astro 05 

Battery Pack 

7 AA Nickel-Cadmium 

Propeller 

Tornado 10-6 



Upon receiving the request for proposals on January 18, 1990 , 

low Reynold's number station keeping flight vehicle, Desig p 

resign phase in pursuit of such a vehicle. Bach-met printed his 

ideas in regards to the mission requirements and goals wh 

crouo in their formal presentation of the Design Requirements and 

EE 

— d tj rr 

propulsion, and structures l presented the status of their 

S £ZL ' — ~ — on March 

28, 1990. 

Upon receiving approval for their design, the group began the manufacturing 

and formal documentation phases. Th. manufachiring began April 17 19 0 

and progressed rather rapidly due to the assigning of a chief 

eagerness of the other member. The aircraf t assmnb^ - »mp.e* ^on Apn. 

29 and was rolled out for ground testing on May 1, 1990. The airc 

29 and was roue » to lhe documentation of the 

flown successfully on May 4, 1990. mg ,. , his fina i 

design each member contributed his expertise and studies to this tin 

written report. It is presented as an overview of this design process 

guide for other design teams to follow. 

The members of Aerospace Design Group D have their 

name to one that is less generic and more suggestmg of their d> ec . 1 

are now referred to as Aerospace Engineers for the wen y 
ir* abbreviated AETC. Likewise, they *-- 5 — 

SCREEM-J4D. This is an acronym derived from t 

member's name: Steve, Chris, Rob, Eric. Ed, Mike, and Jeff. The 4D 

for "For Diploma". It is often referred to as the J4D. 
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This document represents the Final Design Proposal as presented by the 
Aerospace Engineers for the Twenty-First Century in response to the request 
for proposals for a low Reynold's number aircraft. It progresses in sections 
with the first being a look at the preliminary concepts. The next five sections 
deal with the particular technical areas related to the design - aerodynamics, 
stability and control, propulsion, weights analysis, and structures. 
Performance parameters are then considered followed by discussion related to 
the technology demonstrator. Discussion of the design ends with the 
conclusions section and a final word from the design team. The references 
and appendices represent good sources of information for future design 

teams. 


xi 


SECTION I 

PRELIMINARY CONCEPT 
DESIGN 


A. Request for Proposals 


The following Request for Proposals (RFP) provided the design requirements 
for a remotely-piloted vehicle (RPV) to be created by AETC. 


UNIVERSITY OF NOTRE DAME 

DEPARTMENT OF AEROSPACE AND MECHANICAL ENGINEERING 

AE441: Aerospace Design: Request for Proposals 
Spring 1990 

FLIGHT AT VERY LOW REYNOLDS 
NUMBERS - A STATION KEEPING MISSION 

OPPORTUNITY 

Most conventional flight vehicles are designed to operate in a flight regime 
such that the Reynolds number based on mean wing chord are in excess of 10 6 
and some currently approaching 10®. Recently there has been interest 
expressed in vehicles which would operate at much lower Reynolds 
numbers, less than 10 5 . Particular applications are low speed flight at very 
high altitudes, low altitude flight of very small aircraft and flight in other 
planets' atmospheres such as Mars. There are many unique problems 
associated with low speed flight which pose challenges to the aircraft designer 
and which must be addressed in order to understand how to exploit this low 
Reynolds number regime. Since many of the anticipated missions for this 
type of aircraft are unmanned, it is necessary to couple developments in 
unmanned aircraft development with our knowledge of low Reynolds 
number aerodynamics in order to develop an aircraft which can fly as slow as 
possible at sea level conditions. This study will help to better understand the 
problems associated with flight at these very low Reynolds numbers. 
Considering the potential applications, the aircraft must be very robust in its 
control and be highly durable. 
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OBJECTIVES 

1. Develop a proposal for an aircraft and associated flight control system 
which must be able to: 

a. Maintain level controlled flight and fly a closed course at flight 
speeds corresponding to Reynolds numbers less than 2xl0 5 and as close to 
IxlO 5 as possible. The greatest measure of merit is associated with achieving 
the lowest mean chord Reynolds number possible and maximizing the loiter 
time on a closed course. 

b. Be maneuverable and controllable so that it can fly a closed pattern 

and remain within a limited airspace. 

c. Use a propulsion system which is non-airbreathing and does not 

emit any mass (i.e. rocket, etc.). 

d. Be able to be remotely controlled by a pilot with minimal flying 

experience or an autonomous onboard control system. 

e. Carry an instrument package payload which weighs 2.0 oz and is 

2"x2"x2" in size. 

2. Take full advantage of the latest technologies associated with lightweight, 
low cost radio controlled aircraft and unconventional propulsion systems. 

3. All possible considerations must be taken to avoid damage to 
surroundings or personal injury in case of system malfunction. 

4. Develop a flying prototype for the system defined above. The prototype 
must be capable of demonstrating flight worthiness of the basic vehicle and 
flight control system. The prototype will be required to fly a closed figure 8 
course within a highly constrained envelope. A basic test program for the 
prototype must be developed and demonstrated with flight tests. 

5. Evaluate the feasibility of the extension of the aircraft developed under 
this project to high altitude station keeping application for atmospheric 

sampling. 

SYSTEM REQUIREMENTS AND C ONSTRAINTS 
The system designed shall satisfy the following. 

a. All basic operation will be line-of-sight with a fixed ground based pilot, 
although automatic control or other systems can be considered. 

b. The aircraft must be able to take-off from the ground and land on the 

ground. 
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th . i oi t e r time within a restricted 
c The aircraft must . spacing of 150 ft between two 

altitude range on a t gu , 

pylons Which define the course. „ able to be accomphshed 

4 Ground handUng and system operation 

by two people. M be disassembled for transportation 

e The complete aircraft must be ab ^ larger than 2'x2'x4'. 

and storage and fit ^ tWn a " ^ operation are critical. A complete safety 
f Safety considerations for systems op 
assessment for the system is require . 

rt flight tests for the ^ 

Demonstra point on the course. shown on the 

must not exceed 25 any po 15( y takeoff region sho 

c. Takeoff must be accomplished wi 

enclosed figure- needed to complete the 3 comp 

d Loiter time will be based on the 

laps in the air. provisions for estimating altitude an ig 

• “ ' 

speed during the tests. 

based observers. Technology Demonstrator mus 

sr— 

rsr^rr- 7 

i. All FAA and FCC regulations tor p 
must be complied with. 
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B. Mission Evaluation 


In order that the aircraft possess the ability to perform the basic mission, there 
are several additional constraints regarding the aircraft and its mission which 
need to be specified. In order to maintain a lightweight and relative y simp > e 
airframe, no ailerons will be used. All banking and turning will thus e 
accomplished by deflection of the rudder. Due to the possible complications 
of flying into and out of ground effect, it is desirable for the aircraft tobe 
above the influence of ground effect before banking into the first _ 
Therefore, the aircraft must enter and leave all turns at the design cruise 
altitude, losing or gaining as little altitude as possible. This may Require 
elevator deflections and/or throttle setting adjustments. y ' 

aircraft will be flown a, an altitude of 20 feet, which is 5 fee i ^ the 
absolute ceiling of the test area. This gives the pilot a margin (altho g 
small) to make adjustments to save the plane in the event it would 

unpredictedly . 

Combining these two constraints, it is evident that upon " 

should climb quickly to its cruise altitude of 20 feet and remain there fo - the 
duration of the mission. By somewhat arbitrarily selecting the take-o 
distance to be 75 feet (which is one-half of the available runway), it becomes 
necessary for the aircraft to climb at a flight path angle of 15 degrees to reach 
the cruise altitude of 20 feet over a horizontal distance of 75 feet. 

As stated in the RFP, the flight envelope is very constrained, and as such the 
aircraft needs to be very maneuverable. With that in mind, AETC feels ha 
the aircraft must have the capability to complete the turns using a 40 foot 
radius leaving no less than 20 feet between the aircraft and any side boundary 
of the flight envelope during the duration of the mission. Since a large lo. er 

time is the greatest figure of merit, flying the mission a, the ^ 

velocity is a must. But everything is related, and since wing surface area and 
thus weight) go up as velocity goes down, a compromise exists somewhere 
between. It is our goal to find that compromise. 
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C Design Requirements and Objectives 


In order to insure that the minimum goals defined in the RFP as _well 
those specified in the above mission evaluation are satisfied, AETChs 
imposed some additional (or more stringent) constraints on its RPV. They 

are: 

1 Maintain controllability while taxiing. Since the aircraft must be ground 
launched, it is necessary for it to be steered while still on the ground at speeds 
lower than those necessary for the lifting surfaces to prov.de th.sdirection 
control. For this reason, the landing gear must be coupled with the av orncs 
system to avoid the necessity of an additional servo to control the grou 

flcHeve airborne flight after a takeoff run of no more than 75 feet This 
constraint forces the aircraft to be off the ground and climbing 
than half of the runway left. This means that should some 
immediately upon takeoff, there will still be sufficient runwa, 'left to * 
aircraft back on the ground without running it mto a wall. Addit y, 
gives the aircraft a relatively long time to reach its cru.se altitude of feet, 
mandating a flight path angle of approximately 15 degrees. 

3. Climb at no less than 6 feet per second. Based on a cru.se veloa y o fee 
per second; if the aircraft takes off in 75 feet this rate of d.mb will allow the 
aircraft to climb to its cruise altitude of 20 feet by the end of the runway. 

4. Maintain controllable flight at a chord Reynolds number o no more to 
lxlO 5 This is the lower limit of the Reynolds numbers targe m 
Setting the design point at this Reynolds number ensures that should small 
increases in velocity become necessary to avoid stall, the Reyno s num er 
will not be pushed over the 2x1 0 5 upper limit 

5. Maintain altitude during turns of no less than 270 degrees on no more 
than a 40 foot radius at bank angles of up to 30 degrees. Th. con* mm. 
evolved from safety considerations for the operator of the aircraft and for 
preservation of the aircraft itself. The figure "8" course mandates the turns 

at least 270 degrees; without these turns the aircraft will sure y in ‘ sw ^ 
into a wall. The radius requirement ensures that the aircraft will P^s 'he 

ability to turn well within the maximum allowable radius of 60 ee give y 

the geometry of the test flight area. FinaUy, the bank angle constraint ensures 
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that even at severe bank angles, the aircraft will be in no danger of dragging a 
wing on the ground due to lost altitude. Please note that in order to 
maximize the loiter time, the final test flight may not be flown at these target 
values, but the aircraft must possess the ability to fly there should the 
situation or need present itself. 

6. Weigh no more than 3 pounds. This constraint evolved from 
comparisons with other remotely piloted vehicles of similar size or mission. 
The primary consideration here was that a lighter weight would allow for a 
lower velocity and subsequently higher loiter time. 

These constraints provided the guidelines by which AETC's aircraft, the 
SCREEM-J4D, was designed and ultimately constructed. While not 
specifically directed at any one area or sub-system of the aircraft, these 
constraints had far-reaching implications in all aspects of the aircraft's design 
and construction. Table 1-1 is a summary of these constraints. 

TABLE 1-1: Design Constraints for the SCREEM-J4D 


1. Maintain controllability while taxiing. 

2. Take-off in distance <75 ft. 

3. Climb at > 6 ft/s and fly at altitude = 20 ft. 

4. Maintain contollable flight at a chord 
Reynolds number < lxl(A 

5 . Maintain altitude during turns >.270 degrees 

with turn radius < 40 ft at bank angles <30 degrees. 


D. Concept Selection Study 

The selection of a concept became a laborious task for AETC once the Design 
Requirements and Objectives (DR&O) were completed. The constraints given 
by the RFP in addition to those in the DR&O had to be weighed very carefully 
in the process of selecting a final concept. The approach that AETC took was 
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considering three different concepts which were unrelated to each other but 
were proposed by at least one member of the group. This was reasonable due 
to the fact that each concept was considered in a more general form with 
relation to the constraints. A more detailed analysis followed once one of the 
three concepts was selected. 

One of these concepts was a biplane configuration with a conventional 
empennage (Figure 1-1). The main advantage of the biplane was that it could 
provide extra lift which is vital to flying at low Re. Nonetheless, the 
disadvantages were too numerous. The complications in resolving the 
interference effects, the construction and connection of the two wings, and 
the additional weight with the added wing were the main disadvantages 
which prompted our decision to reject the biplane concept. 

Another concept which AETC labored over for about a week was a canard 
configuration with a vertical tail (Figure 1-2). The canard was considered able 
to provide extra lift while maintaining maneuverability and not adding 
weight. However, after doing further research on low speed RPV flight, 
AETC discovered that the canard would have to be too large in order to have 
a significant contribution toward the total lift. As a result, the airplane would 
have just as much if not more wing surface area as a conventional airplane 
without the convenience of the classical stability equations. Thus, the canard 
concept was considered "pretty" but not the most efficient concept for the 
mission. 

A third concept which AETC decided to chose for the mission was the 
"sailplane" type airplane with a conventional empennage. The most 
important factor in the selection of this concept was that the theories and 
equations which were such a big part of undergraduate aerospace study could 
be most readily applied to the conventional airplane. In addition, the 
learning process of the design and construction of this concept is useful in 
gaining a better understanding of the process involved with the 747 s, DC-lCKs, 
etc., becoming a realistic part of the transportation industry. 


Table 1-2 is a comparison matrix of the three concepts with regard to several 
critical factors found in the DR&O or assumed from low speed flight. A +1 
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Figure 1-1: Biplane Configuration 





UJ 

vjj 

a 




denotes that the design would benefit from the concept while a -1 denotes 
that it would be detrimental. This table enabled AETC to be more easily 
convinced that selecting a high aspect ratio conventional airplane as their 
concept would best suit the mission. 


TABLE 1-2 : Comparison Matrix for Concepts Considered 


FACTOR 

BI-PLANE 

CANARD 

CONVENTIONAL 

Added Lift 

+1 

0 

0 

Ease of Construction 

-1 

+1 

+1 

Weight 

-1 

0 

0 

Use of Equations 

-1 

-1 

+1 

Maneuverability 

-1 

-1 

-1 

TOTAL 

-3 

-1 

+1 


The initial analysis for providing a little more detail to the concept was based 
on the Reynold's number constraint of 100,000 and the weight constraint of 3 
pounds. Given that the Re is proportional to the mean chord length and the 
velocity and inversely proportional to the kinematic viscosity, it was seen 
that by selecting a target Re and varying the mean chord, a velocity which 
satisfied the constraint would be produced. Using this velocity and chord and 
estimating a target lift coefficient, the span necessary to produce the required 
lift was calculated. This process was aided by use of Computer Program A 
which is located in Appendix I. However, knowing that no piece of the 
aircraft can be longer that 4 feet required that the span be kept reasonably low 
in order to avoid having a wing with too many sections. Using this method 
with the lift coefficient of 0.9 and a weight of 3 pounds, it was found that the 
allowed velocity and required span for a chord of 5 inches was 40 ft/ s and 4.2 
feet, respectively; while for a 12 inch chord, the same quantities would be 16 
ft/s and 11 ft. 

Because the loiter time is dependent on velocity, the first case is not optimal, 
and because lightweight is of critical importance, the large wing area of the 
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second case, with its smaller velocity, is not optimal. Obviously, AETC s 
solution fell between these cases. 

AETC made seveal other decisions about the conventional concept once it 
became obvious it was best suited for the mission. In addition to the high 
aspect ratio wings, AETC decided against tapering the wing because the chord 
at the tips would become too low and possibly introduce problems with 
extremely low Re flight. Additionally, an untapered wing is much easier to 
construct since the airfoil sections are all the same size. Ailerons were 
discarded as a possibility because the added weight of a servo and the 
construction required outweighed the advantages of better roll performance. 
Past RPV's and research proved that the same roll performance could be 
attained with just the elevator and rudder in the first place. Sweep is not 
advantageous for low speed flight since it reduces the effective velocity e 
wing sees. However, dihedral aids in roll control. Winglets were consider 
in the concept study in order to reduce induced drag. But when the ig 
aspect ratio of the wing was considered, it was decided that the effect o 
winglets would be too small to justify the added weight and construction 
difficulties. A low-mount wing was selected to help destabilize the ro w 1 
would facilitate the figure-8 maneuvers. 

The fuselage was selected to be a box configuration for ease of construction 
and containment of the propulsion system. The aft section of the fuselage 
was tapered to reduce weight and skin friction drag. 

The empennage was designed to contain two control surfaces rudder a 
elevator. A relatively large rudder (55% of vertical tail area) was needed to 
produce the yaw required to attain the desired roll moment. The rudder, in 
conjunction with a large wing dihedral (10 degrees), enabled ‘he aircraft to 
perform the figure-8 turns as desired. The elevator comprised 33% of the 
horizontal tail area, ensuring positive longitudinal stability and control. The 
center of gravity was placed behind the aerodynamic center so that the 
horizontal tail would produce lift rather than a downward force. 

Lastly, the landing gear was comprised of three wheels - one on either sideof 
the fuselage mounted to the wing and one in the rear as a tail dragger. The 
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two on the wing were spaced well enough to prevent tipping during take-off. 
The tail dragger was connected to the rudder so that one servo could be used 
for turning control during ground roll and flight. 

These characteristics of the basic concept may be seen in the three view 
drawing in the Executive Summary. 
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SECTION n 


AERODYNAMIC 

CONFIGURATION 



A. Airfoil Selection 


At low Reynolds numbers, viscous effects become more significant with 
respect to inertial effects. The viscous region extends beyond the surface of 
the body both in the lateral and longitudinal direction. The wake retards the 
oncoming freestream, while interaction between the viscous and inviscid 
regions is strong and nonlinear. Therefore, small changes in the position 
(i.e., attack angle) of the body exposed to the flow will cause significant 
changes in the surface pressure distribution. Due to the effects of this viscous 
region, the flow is less capable of negotiating an increase in pressure. The 
adverse pressure gradient can grow large enough to affect a peculiar 
aerodynamic phenomenon— a separation "bubble," which is an enclosed 
region of stagnant air that develops on the upper surface(s). 


The flow over a smooth airfoil at low Reynolds numbers is often laminar, 
unless means are employed to disturb the flow and thereby induce 
turbulence. (Disturbing the flow energizes the boundary layer, thereby 
increasing inertial effects and delaying separation). As the Reynolds number 
decreases in laminar (undisturbed) flow, a separation bubble develops on the 
upper surface and effectively increases the camber of the airfoil, thus 
providing additional lift. However, the separation bubble is extremely 
sensitive to disturbances (such as surface roughness and sound), and/or an 
increase in the angle of attack. If the separation bubble breaks down, the 
pressure drag will increase drastically. In addition, hysteresis is particularly 
significant in airfoil data taken at low Reynolds numbers and at high angles 
of attack. Therefore, the airfoil performance depends strongly on whether the 
angle of attack is increasing or decreasing. Thus, in selecting an airfoil, it is 
important to keep in mind the characteristics of that airfoil in the low 
Reynolds number regime, as well as the flight plan dictated by the mission 
requirements. In particular, selection criterion should include gradual stall, 
high stall angle and/or low stall speed, and lift performance which suits the 
flight plan and aircraft gross weight. 

Measures of Merit: To select the airfoil intended for operation in a low 
Reynolds number regime, the following characteristics are necessary. 

1. Favorable lift curve slope and high CLmax 
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2. Gradual stall 

3. High lift/drag ratio 

4. Ease of fabrication 

Design Constraints and Consideration,: There are severai fixed aerodynam.c 
performance parameters which are influential in the des.gn of the arrfo.1. 

The following were most critical to the J4D: 

1. Cruise Reynolds number (based on mean wing chord), Re c -iu 

2. Cruise velocity, V oo 23 ft/s 

3. Mean wing chord, c = 0.683 in. 

4. Freestream density, p=2.37 x 10-3 s i U gs/ft3 

p .. H ^ these constraints, there are two are other important guidelines which 
w^e used during the design phase. The first was ^ capability o prec, 
fabrication of the seiected airfoil. For example, very thin mrfods, a g 
popular for low Reynolds number operation, are 
precisely with the materials and tools which were available to AE • 
addition thin airfoil sections are subjected to higher stresses in flight and 
under static conditions. The second important consideration was that st 
characteristics must be amenable to the desired flight “"dit.ons^ Given 
requirements for lift performance, the maximum 1. 1 ^ire^must 
achievable without flying at a geometric angle of attack close 
(a „) Therefore, finite lift characteristics yielded by the airfoil paramete 
n7 he wing geometry must be acceptable for the mission requirements 
while restraining the wing size to reasonable limits to minimize the aircraft 

overall weight. 

a mntiniM* Because the airfoil data used in the selection process [Ref 
" aTfle^ r s nlbers close to, bu, no, exactly, the design operating 

Reynolds null for our vehicle (iO^ interpolation was necessa^ o 

determine the important p—- |£JJ ^timated 

sIodc otcfaii/ and. maximum L/D). Fig . , i j 

data 'for Re= 10?, for the NACA 4415. A linear variation with Re y n0 

number was assumed for determining this data. In addition, hysteresis effe 
wT* eglected since the intended flight plan is a steady state cruise a, a 

constant a. 
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x: z xrsr r — *- - *— 

options. U> pro e tric attack angles well below oxtail- 

delaying separation, and (2) fly g of these options. 

There are, however, disadvantages associated wi 

because the desired lift performance can be achmved angies 
below Ostall. option (2) was employed on the SCREEM-J • 

The research for airfoils at low Reynolds numbers was conducted with 

concern with high stresses in these thin region, tajddrhon,^ ^ 

requirements of laminar ^ disturbances , particularly 

used. (Option V Thus although better lift characteristics 

associated skin friction drag penalty). Thus, g flQW airfoils , 

are produced in the low Reynolds number regi y thicker 

disturbance of the laminar flow is very difficult to avoid; therefore, a 

airfoil profile is desirable. 

After a thorough evaluation of the low ^ ££%£ 

two were selected for final analysis, the AH 79-100C ^ancU ^ ^ ^ 

Ficmres 2-2 - 2-8). A finite wing lifting-line analys P® 

airfoil using the computer program listed as Item 2 ml a[>gle 

computer analysis req ^^“ “ s^plXi* the wing planform 
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Figure 2-2: Airfoils Considered 


Scale Profile ollhe AH-79-100C Airfoil . 
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Figure 2-4. AH 79-100C Finite Wing Slope 
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Figure 2-7. AH 79-100C CDi vs. Geometric Angle of Attack 



Figure 2-8. AH 79-100C Drag Polar 
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lift characteristics and the measures of merit above, it was not chosen as the 
main wing airfoil. The fact that it is so thin introduced construction 
difficulties and posed some concern for the airfoil’s structural stability and 
consistency. Consequently, the AH 79-lOOC's better performance 
characteristics were not convincing enough to accept its disadvantages, and 
the NACA 4415 was selected by default. 

Final Selection: The NACA 4415, based on the information available in 
Reference 12 and interpolated approximations for Re=10 5 , has the following 

approximate parameters at infinite span: 

Lift curve slope: 3.72 per radian 
Zero-lift angle of attack: - 4 degrees 

Qmax : ~ 

L/Dmax : ~ 15.5 

The values for maximum lift coefficient and maximum lift-to-drag ratio did, 
of course, change when a finite aspect ratio was considered. 


B. Wing Geometry 

Important Parameters: To specify the wing geometry, the following 
parameters must be determined: 

1. taper 

2. span 

3. geometric and aerodynamic twist 

4. dihedral (discussed in Section III) 

These parameters were determined using these fixed (design) values: 

1. mean-chord Reynolds number, Re= 10 5 

2. cruise velocity, Voo= 23 ft/s 

3. freestream density, p* 2.37 X 10* 3 slugs/ft 3 

4. freestream viscosity, |i“ 3.76 X10" 7 Ns/m 2 

The wing span and chord variation together determine the planform 
[projected] wing area. The mean chord is fixed at a value of 0.68 ft by the 
cruise velocity and desired Reynolds number. 
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Measures of Merit: The planform geometry was selected based on the 
following characteristics: 

1. finite lift curve slope 

2. stall propagation pattern 

3. spanwise lift distribution 

4. induced drag 

5. structural weight of wing and wing loading 

It is desirable to obtain finite lift characteristics commensurate with the 
performance requirements and overall system weight. To minimize the 
overall system weight, the planform geometry was selected such that the 
overall weight of the main wing was at a minimum, given the performance 
requirements. 

To obtain a favorable stall pattern, the planform geometry was selected such 
that the maximum lift coefficient occurs at the root of the wing. The stall will 
thus propagate from the root. This is particularly desirable for minimizing 
bending moments in the wing. This stall pattern is also desirable if control 
surfaces (i.e., ailerons) are to be located on the outboard wing sections. 
Although the SCREEM-J4D was designed without these control surfaces, 
minimizing wing bending moments is still extremely important. 

An additional means of minimizing the stresses incurred in the wing is to 
obtain a lift distribution which, when coupled with the wing surface area, 
yields a reasonable wing loading. This total wing loading is obtained by 
integrating the spanwise lift distribution. 

Design Constraints and Considerations: As previously mentioned, airfoil 
selection and wing geometry are influenced by several parameters which 
were fixed by AETC in the interest of successfully meeting the mission 
requirements. Table 2-1 is a summary of these fixed parameters. 
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TABLE 2-1: Design Parameters for Airfoil & Wing Geometry 


PARAMETER 

VALUE 

Poo 

2.37 X 1(T 3 slug/ ft 3 

Voo 

23.0 ft/s 

Re c 

It) 3 

aircraft gross weight 

3.0 lbs 


Other constraints and decisions were made during the preliminary design 
concept (in addition to fixed parameters): 

- stall angle - required Cl must be achieved at a < ctstall 

- no taper (taper ratio = 1.0) 

- no aerodynamic twist 

If a taper ratio is employed on a wing intended for operation at low mean- 
chord Reynolds numbers, separation effects will cause substantial lift 
reduction and drag increase at the tip section. This is due to a lower chord 
than the mean chord and hence the local Reynolds number is smaller than 
the mean-chord Reynolds number. Therefore, the taper ratio of the wing will 
remain fixed at 1.0. 

Since the mean-chord is fixed by the Reynolds number and cruise velocity, 
and the chord is constant along the entire span of the wing, the span is 
determined by a choice of the planform area (or vice versa). Obviously, the 
choice of airfoil section also directly affects the finite lift characteristics of the 
wing. The (infinite] airfoil section parameters which affect the finite lift 
characteristics are the lift curve slope and the zero-lift angle of attack. Because 
the lift generated by the wing depends on its wetted surface area, the lift is 
dependent on the choice of wing area or span. In this case, area and twist are 
the only wing geometry parameters which need to be specified to generate a 
certain lift force. Aerodynamic twist will not be employed, in the interest of 
simple wing construction. Thus, our wing geometry options are limited to 
span and geometric twist configuration. As a result, the independent 
variables whose effects on lift have been investigated are the planform wing 
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area, geometric twist configuration, infinite lift curve slope, and infinite zero- 
lift angle. 

An important objective of AETC is to develop a design which is easily 
constructed, yet capable of successfully meeting the mission requirements. 
Thus, the independent variables investigated were chosen in part on the basis 
of their relatively insignificant effects on construction capability. Because 
aerodynamic twist complicates the wing construction, and (as it turns out) the 
desired lift performance can be achieved without employing this technique, it 
has been disregarded as a possibility in the design of the wing. Further, the 
weight of the wing is an important consideration. Minimal weight of the 
aircraft is particularly crucial in low-speed flight. Therefore, the wing size 
must be sufficient to achieve the maximum desired lift force, yet no larger. It 
is important to remember that because interference effects with the fuselage 
are not taken into account in the lifting-line wing analysis, the performance 
predictions of the actual wing will be slightly inferior to those calculated by 
using finite wing theory. 

Assumptions: The analysis that follows is based on several important 

assumptions: 

1. The mean-chord Reynolds number is based on the cruise 
velocity, which is fixed at 23 ft/ s. 

2. The angles of attack used in the analysis are sufficiently small for 
valid applicability of the monoplane [finite wing] equation (Eqn. 

2.58 [Ref. 2], Eqn. 5.51 [Ref. 1]). 

3. Symmetric load distribution. 

4. Effects of fuselage on wing lift performance are not taken into 
account. 

Methods of Analysis and Selection: To determine the span efficiency factor, 
finite lift curve slope, induced drag vs. geometric angle of attack, and 
spanwise lift distribution, a finite wing analysis program was employed at 
geometric attack angles in the range, 1 degree < a < 11 degrees. The finite 
wing analysis program employs the monoplane equation of Prandtl s lifting 
line theory [Ref. 2]. This equation is: 
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a(8) = 


2b 

7C(C(0)) 


sin w0 x 

I (A rt sin n0 + n A„ — — ) + a L =o 


where a(0) = geometric angle of attack 

cxl=o = airfoil zero-lift angle of attack 
b = wing span 

c(0) = functional form of the chord variation (in our case, c(q) is 
constant) 

A„ = Fourier coefficients from the proposed Fourier series 
representation of the circulation distribution. 

To reach a decision regarding the optimum combination of the planform 
geometry parameters considered, the required lift during the take-off and 
cruising phases of the flight were determined based on weight estimation, 
desired rate of climb, and an estimate of the losses incurred due to fuselage- 
wing interference. The desired C L for cruise was found according to: 

W 

C L = T 

ipVoo 2 (0(1) 

where W is the total aircraft weight, and (c) is the mean chord (which is fixed 
by the design). The (finite) lift coefficient required during take-off was also 
computed with this equation replacing W with the required lift during take- 
off. The required lift was determined from the following equation: 


L = Wsin0, 

where 0 is the desired flight trajectory. To account for the interference effects, 
a factor of safety was taken to be 1.2 based on an estimate of the effective loss 
of lift. The values for C L determined by the above equations were multiplied 
by this factor to yield realistic requirements. Note that because the induced 
drag coefficient, C d ,i, (and therefore the total wing drag coefficient) varies 
quadratically with C L , it is desirable to obtain the minimum value of C L 
necessary to provide the required takeoff performance and sustain the aircraft 
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in flight at the desired cruise speed. Thus, a compromise was made between 
lift coefficient and induced drag coefficient. 

Results- Given the measures of merit, the design constraints and the 

RC • • c lu mat a 4415 airfoil (see "Airfoil Selection ), 

performance charactenstics of the NACA 4415 airtoii tsee 

Table 2-2 represents the final wing geometric and aerodynamic paramete . 

TABLE 2-2: Selected Wing Geometry 


Chord, c= 0.683 ft (constant) 
Span, b= 8.0 ft 
Taper ratio, X= 1.0 
Aero. & Geo, twist: none 


Planform area, S= 5.46 ft^ 

Span efficiency factor, e=0.89 
Average wing loading = 0.55 lbf/ft 2 
Finite lift curve slope, a= 0.073/ de 


Analysis>of Figures: The finite lift curve slope for the NACA mS may be 

estimated from Figure 2-3 as 0.0733/degree. It is evident, afto that the 
maximum desired lift coefficient which is equal to 1.2 is "d at an 
incidence angle of -11.5 degrees. The desired cruise lift coefficient, 0.87, 
achieved at -7.5 degrees. As expected, the induced drag coefficient vanes 
parabolically with angle of attack. The spanwise lift distribution for 
NACA 4415 is shown in Figure 2-9 and Figure 2-10 for angles of attackrang ‘ ng 
from 1 to 10 degrees. As is evident from these plots, the dismbutions 

approximately elliptical. 


C. Drag Estimation 

Paramount in any aircraft's design phase is an accurate 

forces produced by the aircraft in flight. Since the power required and many 
other parameters depend almost exclusively upon the aircraft's drag, t is rag 
estimation needs to be constantly updated as the aircraft nears e J*” 
phase The method used here is a combination of one presented in Mr. 
Jensen's thesis "A Drag Prediction Methodology for low Reynold^Numters 

Flight Vehicles" [Ref. 6] and another presented by Dr. Nelson in E 

Fliflhi Mechanics [Ref. 7]. 
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CL (finite span) 


2-9. NACA 4415 Finite Wing Lift Distribution 
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2-10. NACA 4415 Finite Wing Lift Distribution 

(AOA= 6 to 10 degrees) 
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In general, the drag polar for an aircraft is made up of a parasite (or form) drag 
component and a lift-induced drag component. The general formula for the 
drag coefficient is then: 

Cd = Cdo + Cl 2 / eicAR 

where C Do and e are obtained from a method called the drag breakdown and 
AR is the aspect ratio of the wing. 

Determination of C Do : The value for the parasite drag coefficient, C Do , was 
determined by summing contributions from all of the aircraft’s components 
and dividing the result by the reference area of the aircraft, in this case wing 
planform area. In general, this quantity is then boosted by ten percent to 
include interference effects. The general equation used to calculate C Do is 
shown here: 


Cdo = L(Con An)/S re f 

The II's represent the individual components of the aircraft and the Cd and A 
of each II are given by certain parameters pertaining to the shape and position 
of the component. The drag breakdown for the aircraft is summarized in 

Table 2-3. 

TABLE 2-3 : Component Contribution to Parasite Drag Coefficient 


Component 

Con 

An 

(Con An) 

Wing 

.007 

5.46 ft 2 

.03822 

Fuselage 

.110 

.111 ft 2 

.01221 

Vertical Tail 

.008 

.38 ft 2 

.00304 

Horizontal Tail 

.008 

.625 ft 2 

.005 

Landing Gear 

.017 

5.46 ft 2 

.09282 


L 

Sum (Con An) 

= .15129 


From this summation, the preliminary Cdo value (using the above equation) 
comes to .0277, and after adding the 15% interference factor, it grows to .0319. 
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Determination of the Aircraft Efficiency Factor, e: The aircraft efficiency 
factor, e, can also obtained from a breakdown method. This method uses the 
following equation: 


1/e = l/e w + l/eb + 1/eo 

When using this method, e w and eb refer to the efficiency factors of the wing 
and body, respectively, and e D is a "other" fudge factor to account for 
interference effects. Using this method, the aircraft's efficiency factor comes 
to .655. Another method introduced in Reference 6 is based on empirical data 
and is a function of the aspect ratio only. The formula given for this estimate 
is as follows: 


e = 1.78(1 - .045(AR)-68) - .64 


The value calculated from this empirical formula was .713. Since we have no 
basis for judgement between these two factors and both come from reputable 
sources, a decision was made to take the average of these two numbers and 
use the efficiency factor thus gained. Our complete drag polar is then: 

C D = 0319 + C l 2 /.682jcAR 

It is interesting to note here that in the calculation of the parasite drag 
coefficient, the landing gear, because of its large reference area (the same as 
the wing) and relatively high Cq more than doubles the value of this 
coefficient. Therefore, if some effort can be made to reduce the drag of the 
landing gear the aircraft will benefit from a much-reduced drag coefficient. 
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SECTION in 


STABILITY AND CONTROL 


A. Horizontal Tail 


The primary responsibility of the horizontal tail for this aircraft is to provide 
longitudinal static stability while allowing the aircraft to be trimmed at a 
relatively small positive angle of attack. A secondary consideration is that the 
tail must not make the aircraft so statically stable as to be unresponsive to 
control input. This section will discuss the process involved with the sizing 
of the horizontal tail; figures of merit with regard to the primary 
responsibility are the aircraft's trim angle of attack and total lift generated. 
The figures of merit with regard to the second consideration are the stick- 
fixed neutral point and the stick fixed static margin of the aircraft; both of 
which are measures of the static stability of the aircraft. 

Aircraft Moment Coefficients - Trim Angle : The computer program listed as 
Item 3 in the Appendix was created to study the trends involved with the 
longitudinal static stability of the aircraft and as such dealt primarily with the 
moment coefficient of the aircraft about its center of gravity. The program 
was used to compute the contributions of the wing, fuselage, and horizontal 
tail to the aircraft's Cmo and Cm a . Using the method found in Reference 8, 
these components are calculated as shown in Table 3-1 (all quantities needed 
for the following calculations are computed as needed within the program). 

TABLE 3-1: Component Contributions to Cm Q and Cm a 


Component 

Cm 0 Contribution 

Cm a Contribution 

Wing 

Cmow = Cmacw + ^-LwO^cg ~ X^c)/ c 

Gnaw = CLawOQrg - Xac) / c 

Fuselage 

Crrtof = (k2-ki)/36.5ScEwf2(a ow +if)Ax 

Gnof = l/36-5ScZwf 2 (de u /da)Ax 

Tail 

Cmot - *1 VHCLatteo + Hv ~ it) 

Cm a t = -iiVHCLatU - de/da) 


After finding the contributions due to all of these components, the Cmo and 
Cm a contributions were summed resulting in the equation for the moment 
coefficient about the center of gravity of. the aircraft: 

Cm cg = Cmo + Cma a 
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Setting this equation equal to zero will give the angle of attack at which the 
aircraft is "trimmed": 

a = - Cm 0 /Cm a 


Trends & Sensitivities : It should be noted here that the J4D's fuselage 

contributions to the trim angle of attack are constant and almost negligible, 
and although the wing's contribution to the moment coefficient is certainly 
not negligible, it can be treated as a constant for the purpose of the program. 
By assuming that these values do not change significantly as horizontal tail 
parameters change, the program was used to examine the tail's influence on 
the previously mentioned factors. As can be readily seen from the boxed 
equations above, the contributions due to the horizontal tail factor 
considerably into the final output of this portion of the program. Specifically, 
the volume ratio of the tail and the tail incidence angle seem to have a great 
impact on the tail's contribution to the moment coefficients. At this point in 
the design phase; structural, weight, and overall size requirements placed the 
aerodynamic center of the horizontal tail approximately 2 feet behind the 
center of gravity of the aircraft, which is located at 33% of the mean chord of 
the wing. This effectively froze the length to the horizontal tail and made the 
volume ratio of the tail a function of its surface area only. Thus, all of the 
parameters which were previously strong functions of the volume ratio now 
become strong functions of the tail surface area. 

Once the accuracy of the computer program was checked with a series of hand 
calculations, it was modified to calculate the trim angles for a number of 
different tail surface areas and incidence angles. These values generated with 
this program were then plotted to show the variation of trim angle as a 
function of different surface areas and tail incidence angles. As is 
demonstrated in Figure 3-1, the trim angle of attack increases with increasing 
tail area and decreases with increasing tail incidence angle. The multiple 
lines on this plot and the ones that follow are for different values of the tail 
incidence angle. The trim angles shown on the graph are for the fuselage 
reference line of the aircraft. The wing, due to its incidence angle, will see an 
effective angle of attack 4.5 degrees higher than this value. The airfoil 
selected for the wing stalls at approximately 14.5 degrees, so the highest 
possible angle of attack for the aircraft is 10 degrees. The target trim angle of 
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TRIM ANGLE OF ATTACK [deg] 


Figure 3-1 


TRIM ANGLE OF ATTACK CdegD ms. TAIL AREA [ft**23 



TAIL AREA Cft**23 



attack for the aircraft was roughly 6 degrees, allowing close to 4 degrees for in- 
flight maneuvering before encountering problems with stalling the aircraft. 

Horizontal Tail Influence on Lift : The trim angle of attack, though, only 

tells half of the story. Using this angle of attack as weU as the downwash and 
incidence angles, the lift due to the wing and tail were calculated and 
summed. Figure 3-2 shows how the tail area and incidence angle influence 
the lift characteristics of the aircraft. The basic trend is the same: the lift 
increases with increasing tail area and decreases with increasing tail incidence 
angle. Since this analysis is dealing only with the the "steady-level" flight 
phase, the only condition that had to be met here is that the total lift must be 
equal to 3 pounds. Any excess lift is not needed at this point, and any less lift 

will not produce a flyable aircraft. 

Stick Fixed Neutral Point & Static Margin : The second portion of this 

analysis deals with the stick fixed neutral point and static margin of the 
aircraft. The equations for these quantities are given in Reference 8 as: 


Stick Fixed Neutral Point: 

Xnp/c = Xac/c - Cmof/CLaw + nVHCWCWU - de/da) 

Stick Fixed Static Margin: 

SFSM = Xnp/c -Xcg/c 

Although it is difficult to see from the equations standing alone, the variable 
factor having the most influence on these quantities is the horizontal tail 
volume ratio. Knowing that the tail length is effectively frozen, the size of 
the horizontal taU is highlighted as the "most influential factor" in stability 
considerations. 

As a basis for comparison, the horizontal tail surface area was initially sized at 
25% of the main lifting surface area (1.367 ft 2 ) and located 2 ft behind the C.G. 
This combination created a statically stable aircraft capable of producing the 
required lift at a 5.5 degree angle of attack. This configuration, however, also 
produced a Cm a value of -2.6 per radian, along with a neutral point at 100% 
of the mean chord (the end of the wing), giving a static margin of 67% of the 


Stability and Control 


3-3 


TOTAL LIFT GENERATED Clbl 


Figure 3-2 


TOTAL LIFT GENERATED Clb] vs. TAIL AREA 



TAIL AREA Cft**23 


mean chord. While this would certainly be a stable aircraft, it would probably 
prove to be very sluggish. In order to make the aircraft more responsive, it is 
necessary to make the value of Cm a less negative and to move the neutral 
point closer to the center of gravity, thus effectively reducing the static 
margin. Reverting back to the boxed equations, it is easy to see that by 
decreasing the volume ratio (through the tail surface area) Cm a will become 
less and less negative. Figure 3-3 shows the variation of Cm a with the tail 
surface area while Figures 3-4 and 3-5 show how the neutral point in 
percentage of mean chord and static margin are influenced by the same 
quantity. 

Results : With all of these observations in mind, the time for decision 

approaches. The basic criterion remains the same: produce enough lift at a 
relatively low angle of attack. But there was a new twist in the form of the 
Cm a , neutral point, and static margin values which are actual measures of 
the aircraft's longitudinal static stability. Certainly, Cm« must be negative to 
ensure static stability for the aircraft under consideration. Guidelines from 
Reference 8 indicate that static margins of 5% of the mean chord are 
acceptable targets for civilian aircraft. Our preliminary estimate was more 
than thirteen times that value. More constraints were needed to reduce the 
number of choices. Further investigation revealed that the drag increases 
with increasing tail area, and that the lift-to-drag ratio is relatively 
independent of the tail surface area once it rises above 0.8 square feet for all 
tail incidence angles (see Figures 3-6 and 3-7). The overriding factor then 
came from the structures side of the design team, who stated that it would be 
much easier to mount the horizontal tail if the incidence angle were small 
than if it were large. With this information in hand, the decisions were 
made, and the point highlighted on the plots indicate the status of the final 
design. Table 3-2 is a summary of these characteristics of the J4D: 


Stability and Control 


3-4 



ClUalfa 


Figure 3-3 


Cm)alfa vs. TAIL AREA Cft**23 



TAIL AREA Cft**23 
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Figure 3-7 


LIFT TO DRAG RATIO vs. TAIL AREA Cft**23 
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TABLE 3-2: Horizontal Tail Sizing Information 


Volume Ratio 

.356 

Length to Aero Center 

2.005 ft 

Planform Area 

•625 ft 2 

Chord 

.479 ft 

Span 

1.304 ft 

Incidence angle 

-1° 

Cm a 

- .685 rad* 1 

Xnp/c 

.514 

Static Margin 

.184 

Aircraft Lift 

3.01b 

Aircraft a 

5.87° 


As is seen from this table, the static margin is reduced to 0.184, or just over 3 
times the target value of 0.05 set in Reference 8. Additionally, the value of 
C m a has been increased to allow for more control response and less 
sluggishness. The aircraft produces the required lift and is trimmed at an 
angle of attack of 5.87 degrees, allowing just over four degrees for 
maneuvering before stall becomes a problem. Additionally, this tail surface 
area places the aircraft in the region of relatively low drag and high lift to drag 
ratio. 

Elevator : The elevator consists of roughly one third of the horizontal tail's 

surface area. Its surface area is .208 ft 2 and Cm$e is an acceptable -0.683 rad* 1 . 
This elevator area was initially sized based on the rules of thumb found in 
Reference 5. Once the value of the elevator effectiveness was found to be 
"large" enough, this value was accepted and finalized for the design. 

Dynamic Stability : Under its final configuration, the aircraft is dynamically 

stable as well; a damping ratio of .88 ensures steady flight and a time to half- 
amplitude of .105 seconds limits the time the aircraft spends recovering from 
disturbances. The period of oscillation for these oscillations is 1.86 seconds, 
meaning that the motion will damp out after less than one oscillation. 
Clearly this aircraft will be able to perform its station keeping mission in a 
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stable manner using equation found in Reference 8 modeling the aircraft as a 
second order differential equation with constant coefficients. 


B. Vertical Tail 

Lateral stability and control are important considerations in the design of all 
aircraft For this rudder and elevator controlled RPV, the vertical tai an 
wine dihedral were crucial elements influencing its stability, performance 
and handling. Since the SCREEM-J4D was designed without ailerons, the 
rudder was required to induce a sufficient yaw angle, P. such that the vertical 
tail and dihedral would achieve the roll maneuverability usually attributed to 

the ailerons. 


Initial Sizing Analysis: The primary consideration in the sizing of the 

vertical tail and rudder was the ability to induce a side slip angle, or yaw 
angle B (beta). This induced yaw angle, coupled with the design wing 
dihedral, will provide the lateral control for the J4D. The tail surface area 
(S v ), flap effectiveness (x), and the maximum rudder deflection angle (8 r > 
were initially selected comparable to typical RPV values. The initial S v was 
established through the combination of two "rules of thumb:" 


Sv / Sh = 0.33 - 0.40 
S v /S w = 0.07-0.01 


However, values of Sv outside the range of these typical values were also 
considered due to the unusual objective of flying at low Reynolds numbers. 
Table 3-3 contains a further listing of fixed and variable parameters, and the 
typical RPV values used to select the initial values. 
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TABLE 3-3: Vertical Tail Sizing Information 


Parameter 

Value 

Reason 

Surface Area (Sy) 

034 -0.60 ft 2 

RPV typical values 

Position aft CG (l v ) 

1.5 -2.0 ft 

required by horizontal tail 

Rudder Deflection (6r) 

<25° 

stall and structure stability 

Flap Effectiveness (t) 

0.62-0.75 

tail stability, rule of thumb 

Tail Efficiency (n v ) 

1.0 

low wing in propwash 

Lift Curve Slope (Ga v ) 

2k 

elliptical wing assumption 


Lateral Moment Coefficient: The wing-body contribution to Cnp was 

calculated to be -7.21 X 1 O' 5 /deg through an empirical equation and is constant 
throughout the analysis (Eqn. 2.74, [Ref. 8]). This destabilizing contribution 
was negligible compared to the contribution of the vertical tail: 

Cnp v = Vy ny CL(xy (1 + do/ dp) (Eqn. 2.80, [Ref. 8]) 

The value of n v (l+ do/dp) was estimated to be 1.135 and directly 
proportional to S v (Eqn. 2.81, [Ref. 8]). In addition, since the vertical tail is a 
flat plate airfoil, the value of the lift curve slope was assumed to equal 2k. 
The resulting Cnp due to the vertical tail was calculated to be 8.4 X 10" 4 
/ degree, which is roughly ten times the magnitude of Cnp due to the wing- 
fuselage combination. The total yawing moment coefficient was computed 
to equal 7.68 X 10* 4 /deg through: 

Cnp — Cnp v + Cn w f 

The rudder control power, CnSy, was then calculated: 

Cn§ r — -Vy ny x CLqy (Eqn. 2.86, 2.87, [Ref. 8]) 

The ratio of dynamic pressures of the wing and the vertical tail, nv, was 
assumed to remain constant at 1.0, due to the low-wing setting and the 
positive effects of the propwash. From Figure 2.20 [Ref. 8], the flap 
effectiveness t, of 0.72 was selected according the desired ratio of the control 
surface area to the lifting surface area of 0.55. The rudder was sized at 55% of 
the overall vertical tail to afford sufficient control power, without 
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endangering structural stability. The rudder control power, Cngr, was thus 
found to equal -5.33 X 10-4 /degree, at a maximum deflection angle of 25 
degrees. From these parameter ranges, the overall yawing coefficient was 
then determined to vary according to: 

Cn = Cnp |3 + Cn§ r 5r 

In order to solve for the yaw angle, (3, at the trim condition, Cn was set equal 
to zero. Including the trim condition constraint. Beta becomes directly 
proportional to the rudder deflection angle; therefore, the optimum (3 occurs 
with maximum 5r . 

This relation verified and the location of the vertical tail aft of the CG, 
selected in conjunction with the horizontal tail, the vertical tail surface area 
(Sy) was swept through numerous values using the computer program listed 
as Item 4 in the Appendix to achieve a satisfactory design. The vertical tail, 
with tail position ly=1.9 ft and surface area Sy=.3 8 ft 2 (which includes a rudder 
area of .21 ft 2 ), achieved sufficient yawing power to ensure safe completion of 
the mission constraints while remaining in agreement with typical RPVs. 
Numerous values of l v were analyzed and found to have no significant effect 
on p. Since if s maximum effect is about +/- 0.1 degrees, it did not take any 
precedence in the sizing of the vertical tail. The above parameters produced a 
maximum yaw angle of 17.3 degrees at the maximum rudder deflection angle 
of 25 degrees. 

Roll vs. Yaw ; Prior to finalizing the design to the above parameters, an 
analysis relating the yaw to the roll was completed. The computer program 
listed as Item 5 in the Appendix was compiled to compute the roll moment 
due to yaw for various values of dihedral. This code relates the effects of 
wing tip contour and aspect ratio to the roll moment coefficient through: 

Clp = (Clp/Dr + AClp (Table 3.6, [Ref. 8]) 

Since the maximum ordinates of the wing will lie on the mean line, AClp is 
assumed to equal zero. The roll moment coefficient due to yaw angle per 
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degree, (ClR/O, was extrapolated from Figure 3.9 in Reference 8, to equal 
00027. Sweeping the dihedral angles from 4 to 12 degrees, over a range of yaw 
angles of 0 to 30 degrees, a minimum dihedral angle of 7.5 degrees is required, 
as shown in Figure 3-8. Note, the range of required roll moment coe “ ,cie " S 
was initially determined by assuming ailerons. The ailerons assumed in the 
design were selected using "rules of thumb" and RPV's of comparable size. 
Thus an estimate of the minimum required Cl of .035 was computed using: 


Cl = (2 CLctw t da)/Sb J cy dy <Eqn ’ 2 96 ’ * Ref ' ® 

Before ending this analysis, it is important to note that although a decrease in 
tail surface area shows an increase in P. the maximum increase 
approximately +/- 1.0 degree. Still, a smaller surface area does denote le 
weight and less skin friction drag. Thus an optimum design would ‘"elude 
high percentage rudder of a smaller tail area (high T, low Sv). Ta 
summarizes the vertical tail as designed for the J4D. 


TABLE 3-4: Vertical Tail Summary 


Surface Area 
Aspect Ratio 
Height 
Root Chord 
Tip Chord 
Taper Ratio 
Sweep 

Aft CG Position 
Gnp 

Rudder Area 

Max Rudder Deflection 

Flap Effectiveness 


1.38 
8.0 in 
7.2 in 

3.9 in 
0.54 
23.2° 

1.9 ft 

7.68 x 10 4 

0.21 ft 2 
25° 

0.72 

-533 x 10 4 


C. Dihedral 


Four wing 
included a 


dihedral configurations were considered for the J4D. These 
V-shape dihedral and three 3-panel polyhedrals, with various 
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ROLL MOMENT COEFFICIENT Cl 


0.000 


- 0.025 


- 0.050 


- 0.075 


- 0.100 


- 0.125 





panel break locations, as indicated in Table 3-5 and Figure 3-9. Functionally, 
the dihedral links the yaw and roll axes of the aircraft, allowing a yaw 
deflection to cause a rolling moment. This result is due to the increase of 
angle of attack and thus an increase of lift on the forward wing. A comparable 
decrease occurs on the aft wing. The net affect creates a wing generated 
rolling moment. Note however, that the change of angle of attack of the 
forward wing is equal and opposite to that of the aft wing. This equality 
provided the wings with a constant total lift within the unstalled region. 


TABLE 3-5: Wing Dihedrals 


Dihedral Type 

Inboard Angle 

Outboard Angle 

Panel Break 

V-shape 
3 panel 


10.0 deg. 

11.0 deg. 
11.9 deg. 
14.3 deg. 

NA 

25% (1ft) 
37.5% (1.5 ft) 
50% (2 ft) 


Analysis and Results : The change in angle of attack (Aa), for the wing is a 

function of the wing dihedral angle and the yaw input angle: 


Aa = tan- 1 (SIN Yaw x TAN Dihedral) ([Ref. 4]) 

The results of the parametric sweep of the yaw angle from 1 to 25 degrees and 
the dihedral angle from 7 to 21 degrees are plotted in Figure 3-10 (see the 
computer program listed as Item 6 in the Appendix). The change in angle of 
attack was plotted verses yaw angle for various Equivalent Dihedral Angles 
(EDA, discussed below). As demonstrated in Figure 3-10, Aa is proportional 
to both the dihedral angle and the yaw angle. The performance envelope and 
the design EDA are indicated. The envelop and the selected A a are 
constrained by the wing's staU angle, the maximum yaw angle of 17.3 degrees 
and the minimum required dihedral angle of 7.5 degrees (above analysis). 
The required wing angle of incidence for steady-state cruise is 10 degrees. The 
stall angle for the wing at a Reynolds number of 100,000 is 14.5 degrees, 
allowing a maximum Aa of 4.5 degrees. Due to the imperfections incurred 
during manufacture, a maximum Aa of 4.0 degrees was selected. 
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DELTA ALPHA (degrees) 


Figure 3-10 


ANGLE OF ATTACK vs. YAW ANGLE for various WING DIHEDRAL (deg) 



YAW ANGLE (degrees) 


O O' 


Given the yaw angle of 17.3 degrees, corresponding to the maximum rudder 
deflection of 25 degrees, an EDA of 10 degrees was selected from Figure 3-10. 
EDA's of up to 15 degrees were also considered for this design, since they too, 
satisfied the envelope constraints. However, the larger EDA’s, while slightly 
increasing the rolling rate, would decrease the wing's lift efficiency slightly 
and increase the structural instability. The EDA of 10 degrees was also 
ascetically more pleasing, thus increasing the J4D's marketability. While it 
would be optimum to reduce the yaw angle, reducing the drag on the vertical 
stabilizer during maneuvers, with a greater design EDA, 10 degrees was 
selected to ensure an unstalled wing. Additionally, while the performance 
requirements necessary for this mission do not "push" the limits of the 
performance envelop for an EDA of 10 degrees; the V-shape dihedral 
configuration was selected due to the minimum outboard angle, to reduce the 
chance of tip stall, should the performance requirements be increased. 

Roll Moments and Roll Rates : The maximum Aa is shown to be 3 degrees 

at the maximum yaw angle of 17.3 degrees for a 10 degree EDA. The change 
in lift over the individual wing was computed, based on a constant lift curve 
slope of 0.1 /deg within the wing's unstalled region. The total moment was 
then summed from the section roiling moments computed from the change 
in lift over 12 panels for an individual wing. Note, due to this dependence 
on the change in lift, the rolling moment of the dihedral wing is 
demonstrated to be independent of the overall angle of attack. A parametric 
sweep of the rolling moment verses the total angle of attack for the leading 
wing was computed (see the computer program listed as Item 7 in the 
Appendix). The roll moment generated by the maximum Aa of 3.2 degrees 
is shown to equal 44.8 ft lb. 

The rolling moment was then non-dimensionalized by computing the 
moment fractions at their semi-span locations. Holding the outboard angle 
(ODA), constant at 10 degrees, for reasons stated above, the rolling moment 
percentages achievable at a maximum Aa were computed for the various 
configurations (see Table 3-6). These rolling moment percentages correspond 
directly to the EDA’s of the various configurations according to Reference 4. 
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Rolling Moment (%) x Outboard Dihedral Angle (degrees) = EDA (degrees) 


TABLE 3-6: EDA vs. Roll Moment 



Moment Fraction 

EDA 

Roll Moment 

V-shape 

1.0 

10.0 deg 

44.8 ft lb 

3 panel 

.91 

9.1 deg 

40.8 ft lb 


.84 

8.4 deg 

37.6 ft lb 


.70 

7.0 deg 

31.4 ft lb 


Note, the ODA is equal to the EDA for a V-shaped dihedral, since the panel 
break occurs at centerline. The rolling moment for the V-shaped dihedral is 
produced over the entire wing, while the flat center panel of the 3-panel 
configurations makes no contribution to the roll power, tasking the outboard 
angles to generate the lift required to roll. Thus, for a constrained wing 
loading and outboard dihedral angle, the V-shape dihedral, with the largest 
EDA, will allow the greatest roll power. This is further demonstrated in Table 
3-7, where given a constant yaw angle, the wing with the greatest EDA 
develops the greatest roll rate. 

TABLE 3-7: ED A vs. Yaw Angle 


1 Yaw Angle 17.3° 

EDA 10° 


EDA 

Roll Rate 

Yaw Angle 

Roll Rate 

10 deg 

28.8 deg/s 

5 deg 

7.4 deg/s 

9.1 deg 

26.2 deg/s 

10 deg 

16.2 deg/ s 

8.4 deg 

24.2 deg/s 

15 deg 

25.0 deg/s 

7.0 deg 

20.2 deg/s 

20 deg 

31.7 deg/s 


The roll rate performance was further analyzed with respect to yaw angle for 
the design wing dihedral of 10 degrees. The data from this analysis is also 
included in Table 3-7. 
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Turning Radius : This coupled relationship was further considered for 

steady-state turning. Constrained to a 40 ft turning radius at the maximum 
bank angle of 30 degrees, the aircraft was found to require a yaw angle of 5 
degrees. Note, this is the required yaw angle to maintain a steady-state turn at 
a turning radius of 40 ft. By increasing the yaw angle to its maximum of 17.3 
degrees, the turning radius of the J4D can be minimized to 28.6 ft. 
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SECTION IV 


PROPULSION SYSTEM 



A. Motor Selection 


The first task in sizing a propulsion system was to select a motor. Since 
design constraints prohibited the emission of particles into the atmosphere, 
liquid-fueled engines of the piston or ducted-fan type were eliminated 
possibilities. Attention then turned to a propulsion system featuring a 
propeller powered by an electric motor, which met the emissions constraint. 
The task then was to select an electric motor. 

The consideration of electric motors involved the Astro 035, Astro 05, Astro 
15, and several generic "others". By generic "others" it is meant those motors 
for which only the maximum power output, weight, size and RPM range 
were known. No information regarding their operating characteristics were 
available. These motors were typically discovered in modelling magazines or 
manufacturers' literature. The motor selection process to be discussed is 
summarized in Table 4-1. The first criterion upon which the selection was 
based was power production. From initial estimates of motor power 
requirements it was realized that the Astro 035 was inadequate. Its power 
production was too low and hence eliminated from consideration. The 
generic "other" engines all appeared attractive in terms of producing enough 
power to meet requirements. However, since no information as to their 
operating characteristics was available, any analysis and performance 
predictions with them would be extremely difficult. This problem could have 
been resolved through inquiries with the respective manufacturers or tests in 
a laboratory. Time invested in either of these practices could have been cost y 
and possibly unrewarding. Due to time and practicality considerations, those 
motors for which the operating characteristics were unavailable were 

eliminated from consideration. 

Attention then focused on the Astro 05 and Astro 15, whose operating 
characteristics were known and available in tabular form. This allowed for a 
performance analysis once a motor had been selected. In regards to 
selection, both produced an adequate amount of power; in fact, well above 
that which initial estimates required. The size of both motors was attractive^ 
That is, the dimensions were such that both were compatible with the size o 
the fuselage and would mount easily to the front of it. The deciding factor 
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became the motor weight. The Astro 05 weighed an ounce less and its 
standard (recommended) battery pack weighed 9 ounces less. Using smaller 
battery packs than those recommended because of lower power requirements 
would lessen the system weight for both. However, the Astro 05 would still 
weigh less and therefore it was selected as the motor to power the aircraft. 


TABLE 4-1 : Motor Selection Considerations 



Astro 035 

Astro 05 

Astro 15 

"Others" 

Power Production : 

low 

adequate 

adequate 

adequate 

Specifications : 

NOT 

available 

available 

available 

NOT 

available 

Size : 

good 

good 

good 

-- 

Weight : 

(Motor /System) 

4.5/ 11.0 oz 

65/ 16.0 oz 

75/ 25.0 oz 

— 


B. Propeller Data Generation 

In order to properly analyze the performance characteristics of different 
propellers, it was found necessary to generate propeller performance 
characteristics (thrust and power coefficients, Ct, Cp) as a function of advance 
ratio J. Since performance estimates on hobby propellers were not available 
and time did not permit actual wind tunnel testing, a data base was created 
using a program entitled "Propeller Performance Analysis for Small 
Computers," written by Barry N. Young [Ref. 11]. This code enabled one to 
generate theoretical propeller performance characteristics based upon simple 
blade element theory. Input to the program were airfoil data (of the propeller 
blade section), estimated flight conditions, and the propeller geometry. These 
geometrical quantities were the blade angle, chord, and thickness measured at 
nine radial positions. 
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Due to the lack of propeller data, it was assumed that all of the propellers used 
a CLARK Y airfoil section which is used for the TopFlight family of hobby 
propellers. It was also assumed that all of the propellers studied had a 
constant blade angle distribution. Which dictated that the only geometrical 
difference between propellers of the same diameter but different pitch was the 
blade angle at each radial position. Thus, for a family of propellers with the 
same diameter, the measurements only had to be performed once, with the 
blade angle (|3) determined by the pitch (P) according to the relation 

p ■ tajl " db 

where r is the radial position measured from the hub (P and r have the same 
length units). 

The propellers studied ranged in diameter (8, 9, 10, 11, 12, and 13 inches) and 
pitch (4, 6, and 8 inches). Each of these diameters were studied with all three 
pitches, hence, the total number of diameter-pitch combinations resulted in 
the analysis of 18 different propellers. A sample of the input geometric 
characteristics and parameters of Barry Young's program, in addition to a 
sample of the theoretical performance results are listed as Item 8 in the 
Appendix. 

Simple blade element theory (utilized in this computer program) helped 
create a data base of performance characteristics for several propellers, but it 
did not give any indication as to how these theoretical results compared to 
experimental results. A comparison was made using experimentally derived 
performance data available for a 13-7, 14-4, and 14-8 propeller. The 
geometrical parameters were measured for these propellers and input into 
the performance program to obtain theoretical results. A study of the 
experimental and theoretical results for the same propellers yielded a definite 
trend. In general, the theoretical results for the power coefficient were much 
lower than the "real" experimental' results. On the other hand, the 
theoretical results for the thrust coefficient were much higher than the 
experimental results. In order to correct the theoretical data, a quantification 
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of these differences was needed. Correctton factors 

ratio (I) were ascertained for both the power coefficient (Cp) d S 

coefficient (Cq). This correction function incorporated the average di 

^wln the ^theoretical and experiments! resuits for bod. .he ^thrus, and 

power coefficients over a range of advance ratios from 0.2 to 

correction function was then added to the theoretical performance 
correction Figures 4-1 and 4-2 compare the 

rharacteristics to yield near real data, rigures 

theoretical performance curves with the corrected performance ™ rves 
" peiur. Tabie 4-2 gives the form of the correction function and the 

appropriate coefficients. 


TABLE 4-2: Correction Functions for Propeller Performance 


Functional Form : 

cf = co + crj 1 + C2*J 2 

+ C3*J 3 

Thrust, Ct '• 

CO = -0.035761 
Cl = 0.085253 

C2= -0.082838 
C3= -0.039541 

Power, Cp : 

CO = 0.022765 

Cl = -0.14785 

C2= 0.32328 

C3 = -0.2654 

Efficiency, t\ ■ 

CO = 0.83487 

Cl * -7.4942 

C2= 19.653 
C3= -17.739 

CT is added to theoretical value; J 

is advance ratio 


in correcting the theoretical data, it was assumed that the correction factor was 

sr£i:sr: rs ? 

performance curves were augmented by these correction unc 10 
greedy enhanced the accuracy of the propeller selection process. 


C. Battery / Propeller Selection 


The next task in sizing the propulsion system was » 
battery pack and propeller to be used with the As o 
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Figure 4-1 

Thrust coefficient vs. Advance Ratio (10-6) 



J 


* — Theory 

0 — Theory with corrections 

Figure 4-2 

Power coefficient vs. Advance Ratio (10-6) 




propulsion system components were not independent of each other, rather 
significantly related. Therefore, the process became one of considering 
numerous battery pack/propeller combinations. Considere were a ery 
packs ranging from four to ten 1.2 volt Nickel-Cadmium batteries, 
battery type would be determined later from a consideration of mill.amp- 
hours needed for endurance purposes. 

In selecting the battery pack/propeller combination the three crucia 1 phases of 
flight - namely, take-off, climb, and cruise - were examined. Each thes 
flight phases were analyzed separately in order to select the battery/prope.ier 
combination to yield the best overall performance. 

Several constraints governed the selection of the propulsion system. These 
are summarized in Table 4-3 and are discussed —out this sec bon m 
regards to the different components and the processes by y 

selected. Firstly, the maximum distance allowed for groun ro on a 
was 75 feet and was set by the Design Requirements and Ob|ectives s Y- 
maximum current allowed in the circuitry was to be 20 — . ^ arose 
from the fact that the Astro 05 motor would be damaged i 
current exceeded this value and it was “enforced” by inserting a 20-amp use 
into the circuitry so that it would blow and safe the motor, 
obviously undesirable to choose a battery 

would draw a current exceeding 20 amps at any time during flight because 
Z would cause an immediate loss of all power. Another constraint was 
minimum power available output of 37 Watts. This would ensure 
plane climb at the desired rate of climb. 

TABLE 4-3 : Constraints Governing Propulsion System 


Maximum Take-Off Distance 
Maximum Current Draw 
Minimum Power Available 
Maximum Battery Pack Weight 
Maximum Endurance 
Maximum Range 


75 feet 
20 amps 
37 Watts 
11 ounces 
78 seconds 
1776 feet 
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The selection of the battery pack and propeller was accomplished in two 
related manners. One method examined the take-off phase and the associated 
constraints. This method examined a wide range of battery pack /propeller 
combinations and narrowed it down to one which was acceptable. The 
second method examined the climb phase. It considered fewer propellers and 
by applying relative constraints eliminated all but two combinations. 
Considering results of both methods, the battery pack /propeller combination 
was able to be selected with confidence. 

Take-off Phase Method : This analysis proceeded with the take-off 

performance of the airplane. This phase was aided by the computer program 
"Take-Off Performance," written by Dr. S. M. Batill [Ref. 3]. It incorporated 
the previously developed thrust and power coefficients as a function of 
advance ratio, the lift and drag characteristics of the airplane, the wing area 
and airplane weight, the ground friction coefficient of the mission flight 
center, battery voltage and resistance, and the torque and current 
characteristics for the Astro 05 motor. With these inputs, the program 
integrated the ground roll performance at 0.05 second intervals and yielded 
the following output: the distance required to achieve the flight speed; 
maximum current draw; maximum power required to turn the propeller; 
time to achieve take-off; and battery power consumed in ground roll. The 
final selection was based upon take-off distance and maximum current draw, 
even though we desired that the maximum battery drain not to exceed 5% of 
the total milliamp-hour capacity. A sample of an input and output data file 
are presented as Item 9 in the Appendix. 

Using the propeller data, each propeller was analyzed for take-off 
performance in conjunction with 5, 6, 7, and 8 batteries, with each battery 
representing 1.2 volts. Although the amount of batteries used would have an 
effect upon the battery pack weight, the weight of the airplane was assumed 
to remain constant at 3 pounds. This weight estimation is based on 7 
batteries, and hence if the number of batteries were less, resulting in a 
decrease of total airplane weight, the overall take-off performance would 
improve. This assumption was validated through use of the program and 
initial hand calculations. 


Propulsion System 


4-6 



Figures 4-3 through 4-5 represent the results determined through use of the 
'Take-Off Performance" code. After analyzing the data, it was found that the 
take-off distance and the maximum current draw defined a definite trend. 
The grid upon which the curves are drawn represents all combinations of 
propeller pitch and diameter studied. Each intersection point indicates one 
such combination with the abscissa indicating the propeller diameter while 
the ordinate is the propeller pitch. As indicated on the grid, several 
combinations were eliminated because they did not overcome friction, take- 
off in under 300 feet, achieve take-off velocity, or produce enough power. 
These are indicated by the asterisks (*) according to the key at the bottom of 
the figure. 

The take-off distances formed parabolic curves that opened to the right 
(Figure 4-3), while the maximum current draw formed parabolic curves that 
opened to the left (Figure 4-4). This trend in data existed for all of the battery 
pack sizes that were studied. In addition to these results, the battery drain in. 
mahrs is represented for the different propeller sizes (Figure 4-5). The data 
does not form any analytical trend because the battery drain is a function of 
current draw as well as total time for take-off. 

Referring to Figure 4-3, it is shown that for 7 batteries, there existed the 
greatest number of propellers which met the take-off distance constraint. 
These being the 10-6, 10-8, 11-6, 11-8, 12-4, and 12-6, all of which yielded a take- 
off distance of 50 feet or less. Figure 4-4, however, shows that only three of 
these propellers fall below the maximum allowable current draw, these being 
the 10-6, 10-8 and 11-6 drawing approximately 15 amps, 19 amps, and 19 amps 
respectively. The 10-8 and 11-6 were rejected because they were considered to 
be within the limit of uncertainty associated with the accuracy of this analysis. 
Hence, this left only the 10-6 propeller to be examined. The 10-6 yielded a 
take-off distance of approximately 45 feet which fell well below the maximum 
allowable take-off distance of 75 feet. In addition to having a relatively short 
take-off distance, this battery pack /propeller combination drew a maximum 
of only 15 amps, 25% less than the maximum allowable current draw. It 
should also be noted that this combination, as shown by Figure 4-5, only 
consumes approximately 15 mahrs during take-off, which is quite affordable. 
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Figure 4-3 

Take Off Distane (ft) 
7 batteries / 8.4 volts 



• » 


Friction Exceeds Static Thrust 
Max Power Exceeds Max Allowable 






Figure 4-5 


Battery Drain (mahrs) 
7 batteries / 8.4 volts 



* Friction Exceeds Static Thrust 
** Max Power Exceeds Max Allowable 





In summary, the take-off performance analysis indicated that the best 
propeller /battery combination suited for take-off, allowing the airplane to 
take-off in under 75 feet and drawing less than 20 amps, consisted of the 10-6 
propeller. 

Climb Phase Method : Initial calculations showed that the climb phase 

proved to be the most demanding in terms of power requirements. From 
estimations of cruise power requirements and desired rate of climb, we were 
able to determine the amount of available power needed to be supplied to the 
aircraft in order to carry out a successful climb. This needed power available 
became one constraint which the battery pack/ propeller had to fulfill. 
Another constraint came from the allowable armature current. As 
mentioned earlier, the Astro 05 motor would be damaged if it experienced a 
current of greater than 20 amps. For the climb and cruise analysis, an 
Airscrew 9-6, a Tornado 10-6 and a TopFlight 12-6 were considered so as to 
represent the desired range of propellers from which to choose. 

A plot of power available against the number of batteries for the three choices 
of propellers is shown in Figure 4-6. From this figure, several combinations 
are immediately ruled out since they do not provide the 37 watts of available 
power needed in climb. The Airscrew 9-6 propeller with a battery pack of 7 or 
fewer batteries fails to produce enough power. Likewise, the Tornado 10-6 
with 6 or fewer and the TopFlight 12-6 with 5 or fewer batteries do not meet 
the power constraint. These combinations were therefore eliminated as 
possibilities. The armature current constraint served to exclude several other 
possibilities. Figure 4-7, which shows the relationship between armature 
current and the number of batteries for the propellers under consideration, 
illustrates these excluded combinations. The Airscrew 9-6 propeller with a 
battery pack of 10 batteries violates the current constraint. The combinations 
involving the Tornado 10-6 with 10 and 9 batteries exceeded the allowable 
current and were ruled out. This same propeller with 8 batteries was also 
ruled out since the current draw of 19 amps was, as mentioned earlier, within 
the limit of uncertainty associated with the accuracy of this analysis. The 
combinations of the TopFlight 12-6 propeller with 10 to 6 batteries were also 
eliminated from consideration due to the excessive current draw. 
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Figure 4-6 

Power Available Variation with Batteries 
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Figure 4-7 

Armature Current Variation with Batteries 
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Summarizing the results of these constraints, only two battery pack/propeUer 
combinations satisfied the requirements. These were the Airscrew 9-6 wi th 8 
batteries and the Tornado 10-6 with 7 batteries. The former was 
through a take-off analysis since it was unable to generate enough thrust to 
exceed static friction. This enabled selection the Tornado 10-6 propeller and 
battery pack of 7 batteries to be used along with the Astro 05 motor. 

The propulsion system had then been specified with the exception of two 
JaraLers. These' were the battery type (AAA. AA. C, or D) and -he mann 
in which the batteries were arranged (series or parallel). The factors 
influencing the battery type were capacity (in terms of m.ll.amp-hourslan 
weight. This discussion will be presented later. The reasons for winch the 
batteries were arranged in series wiU now be discussed. Smce the powe 
produced by the motor is directly proportional to the open circu, voltage 
was advantageous to arrange the batteries as to achieve a high vo ttage. Jlus 
was accomplished by a series configuration since voltage a 
Secondly, the current does not add in series as it would in parallel. There , 
a series arrangement would keep the current low which is a desirable fea ur 
For these two reasons, the batteries were arranged in senes. 


Considering the Tornado 10-6 propeller with 7 batteries. Figure - s 
the current draw during climb at full throttle is approximately 16 amps. 

4 second climb, the power consumed is approximately 20 ma , w l is qu 
affordable. Upon reaching a design altitude of 20 feet, the motor was 
back (by adjusting the voltage with a speed controller) so .ha, the power 
available would be equivalent to the power required. e P ow " . 

referred to is that which corresponds to a steady-level turn with a bankangl 
of 30°. This is so that the airplane would not lose altitude in the turn, igu 
4-8 shows the relationship between armature current and power available fo 
a range of throttle settings. The power required for a steady-level turn was 
about 8 watts. This figure shows that the corresponding current draw 
roughly 6 amps. Although the power required in steady, wings-level flight is 
slightly less than that required for a steady-level turn, the two powers wi e 
considered equal so as to simplify the calculations and produce a conservative 
estimate. Considering the power consumed in takeoff (15 mahr), the power 
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Figure 4-8 


Relation Between Armature Current and Power Available 
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consumed in climb (20 mahr), the current draw in cruise (6 amps), and the 
aforementioned simplification, the endurance of the aircraft is 221 seconds for 
a battery pack of AA batteries and 129 seconds for 2/3AA batteries. The AA 
batteries were chosen because of the large cushion in capacity. Although the 
2/3AA batteries provided enough power capacity for the anticipated mission, 
there were still unknowns that had to be considered. The larger capacity AA 
batteries were chosen because it provided a margin which would allow for an 
aborted take-off, ground roll testing prior to flight, and a safe margin to 
account for standard unknowns. This comfortable increase in capacity with 
the AA batteries also meant an increase in weight, from 2.94 to 5.46 oz. 
Weight estimations for the initial design included a battery pack weight of 9- 
11 ounces. All other parameters in the aircraft design were able to be 
achieved assuming this battery pack weight. Therefore, by using the seven 
AA batteries the weight was already well below initial estimations. Though 
the additional weight savings with the 2/3AA batteries would have been 
attractive, it was regarded as unnecessary in light of this already "low" battery 
pack weight achieved with the AA batteries. The small weight penalty was 
willing to be accepted in light of the increased endurance. 


D. Speed Control 

As mentioned previously, the propulsion system incorporated a speed 
controller in order to throttle back the motor for the cruise portion of flight. 
The take-off and climb phases of flight were performed at full throttle to 
ensure the desired performance. Throttling back for cruise allowed the pilot 
to maintain steady level flight through the turns as well as reducing battery 
consumption in order to increase endurance. Due to the small margin for 
error in throttling back for cruise, it was desired to obtain some idea of the 
degree of throttling necessary so that the pilot might be informed. Figure 4-9 
shows the relationship between throttle settings and armature current for the 
chosen Tornado 10-6 propeller. As revealed by the figure, for a cruise 
armature current of 6 amps (with 7 batteries) the throttle setting was to be 
approximately 50%. With this information, the pilot knew, upon reaching 
cruise altitude, to throttle the speed controller to 50% power. From there. 
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only minor adjustments were required to maintain a steady-level turn and 
then steady-level flight. 
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SECTION V 


STRUCTURES 


A. Structural Configuration 


Before estimation of the center of gravity could be completed, the structural 
configuration of each component was determined- The structural design was 
based on existing plans of other flight worthy aircraft found in modelling 
magazines and other such literature. Regular design status reports by the 
other design groups and aircraft design expert Mr. Joe Mergen also proved 
quite helpful in determining the structural layout for the SCREEM.-J4D. 

Fuselage : The fuselage is mainly a truss structure with minimal planking 

to keep the weight low and to provide added strength. The structure is 
shown in Figure 5-1. It is comprised of primarily balsa wood with birch 
plywood used as mounting plates. The fuselage is 37 in length and has a 
maximum cross section of 4"x 4". The main section of the fuselage was 
designed in such a fashion as to allow for flexibility in placing the internal 
components. This was done so that it would be easy to adjust the 
components if a center of gravity shift was needed. 

Wing : The wing is pictured in Figure 5-2. It consists of 3 sections: one 36" 

inboard and two 30” outboard sections. The wing was designed in this way 
simply for storage purposes (within a 2'x2'x4’ compartment). The design is 
fairly conventional. False ribs located at the leading edge are used for added 
support for the skin. The analysis for the spar configuration is discussed in 
the next section. The main features of the design are the connection of the 
inboard and outboard wing, and the mounting location of the wing to the 
fuselage. Initially, a 'micro-bubble /epoxy' matrix with an aluminum or brass 
rod set in the spar box was considered for this purpose. But a simpler spruce 
connection was decided upon for lighter weight and ease of construction. 4 
spruce beams will be set in the 'box' formed by the spars and webs (see Figure 
5-3), one in the inboard section and one in the outboard section. The 
connection will be secured together with tape wrapped around the wing. 
Mounting of the wing onto the fuselage will be accomplished using the bolt 
technique. Bolt locations are at the forward and trailing edge of the wing 
center. Added support will be provided by securing the wing with a rubber 
band, if needed. 
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Figure 5-3 


Wing Connection 




Empennage : The empennage (vertical and horizontal tails) are pictured in 

Figure 5-4. The design is truss-like to provide strength and low weight. The 
empennage will be permanently fastened to the fuselage. 


B. Spar Configuration 

Before analysis of the spar was performed, the spanwise distribution of shear 
and bending moment for the cruise flight condition was determined. This 
was done by means of the computer program listed as Item 10 in the 
Appendix. Input for the computer code included: 

a. Wing planform geometry 

b. Aerodynamic section properties as a tabular function of spanwise 
position 

c Structural weight per unit span as a tabular function of spanwise 
position 

d. Aircraft dynamic pressure and angle of attack 
which are read from a data file. A second computer program listed as Item 11 
in the Appendix was utilized to determine the centroid location, section 
properties relative to the centroid, and the direct stress due to both bending 
and axial loads for the idealized airfoil section shown in Figure 5-5. Each 
component of the airfoil was idealized as a "lumped area". The maximum 
direct stress, O = ±77 9 psi, was located in the main spars (areas 3 and 8) which 
are spruce. The maximum allowable stress for spruce is 6200 psi which shows 
that the main spars are safe from failure. Even when considering a factor of 
safety of 2.0, the maximum direct stress is <J = 1560 psi, which is still safe from 
failure. Table 5-1 summarizes the section properties, loading conditions, 
direct and allowable stress for each "lumped area". 
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Figure 5-4 
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I Empennage Structural Configuration 



Vertical Tail 



Horizontal Tail 
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Figure 5-5 


Spar Configuration 
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Figure 5-6 


Spar Configuration 
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TABLE 5-1 : Direct Stress Distribution at Wing Root 


Ycentroid = 

0.00 in 

Iyy = 36.59 in^ 

P = 0.00 lb 

^centroid * 

-2.81 in 

Izz = 0.42 in 4 

My * 2.42 in-lb 



I yz = 034 in 4 

M z = 30.24 in-lb 

i 

CTxx(i) psi 

CTxx(all) psi 

Material 

i 

-17.28 

±400.00 

Balsa 

2 

0.09 

±400.00 

Balsa 

3 

-779.12 

±6200.00 

Spruce 

4 

-0.09 

±400.00 

Balsa 

5 

-36.62 

±400.00 

Balsa 

6 

-1.96 

±400.00 

Balsa 

7 

18.37 

±400.00 

Balsa 

8 

779.12 

±6200.00 

Spruce 


An analysis of the joints connecting the inboard and outboard wings was also 
performed. Figure 5-6 shows the configuration of the spars and Table 5-2 
summarizes the data. As the table indicates, the maximum allowable stress is 
not reached in any of the members. 


TABLE 5-2: Direct Stress Distribution at Y=18 in 


^centroid 

= 0.12 in 

l yy 

= 13.14 in 4 

P = 0.00 lb 

^centroid 

= -0.99 in 

Izz 

= 0.44 in 4 

My = 0.85 in-lb 



Iyz 

= 0.11 in 4 

M z = 10.60 in-lb 

i 

a xx (i)psi 


a X x(aU)psi 

Material 

1 

-258.07 


±6200.00 

Spruce 

2 



+400.00 

Balsa 

3 

-7.64 


±400.00 

Balsa 

4 



±6200.00 

Spruce 

5 



±6200.00 

Balsa 

6 

-61.78 


±6200.00 

Balsa 
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C Materials Selection 


The critical characteristics of the materials under consideration were: 

a. Strength 

b. Weight 
c Cost 

d. Machinability 

e. Availability 

Wood, metal, composites, and ceramics were considered and the relative 
rankings of these regarding the critical characteristics are summarized in 
Table 5-3. 


TABLE 5-3 : Relative Material Characteristics 


Material 

Strength 

Cost 

Weight 

Machinabilitv 

Availability 

Wood 

4 

1 

2 

i 

i 

Metal 

2 

2 

4 

2 

2 

Composites 

1 

4 

1 

3 

3 

Ceramics 

3 

3 

3 

4 

4 


The table shows that wood is the best material to be used. Though wood 
ranked last in strength, it still would provided sufficient strength to meet the 
stress requirements mentioned earlier. The bulk of the structure was 
comprised of balsa wood with spruce and birch plywood used where extra 
strength was needed. The skin of the J4D was "Micafilm" manufactured by 
Coverlite. It is a non-woven fiber that "offers low weight, resistance to 
tearing, and good modulus" fit for models where weight control is a primary 
concern. 


D. V-n Diagram 

The V-n diagram is shown in Figure 5-7. First notice that the gust loads do 
not affect this diagram because of the indoor operating environment which is 
free from any significant gusts. There are no negative load factors because the 
aircraft is not designed to operate at a negative wing incidence. The 
maximum operating velocity, V max , of the aircraft is 58 ft/s, the minimum 
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(stall) velocity, V m i n is 18.9 ft/s, and the cruise velocity, V cruise , is 23 ft/s. The 
limit load factor, ni im , is 1.15 (from a maximum bank angle of 30 degrees). 
The positive stalling portion of the curve is found using standard sea-level 
conditions (p = 0.00237 lb /ft 3 ), a maximum lift coefficient of 1.3 and a wing 
loading of 0.55 lb/ft 2 . The ultimate load factor, n uU = 1.73, is based on a factor 
of safety of 1.5 (n uit = 1.5n, im ). Special attention should be made to the fact 
that the cruise point (n=1.0, V=V cru ise)/ pushes the limits of the flight 
envelope. Therefore, the aircraft must not be flown at an attitude such that 
disastrous loads are placed on the aircraft. The factor of safety, however, does 
allow for somewhat significant loads to be placed on the aircraft, if necessary. 

TABLE 5-4 : Special Flight Conditions Illustrated in V-n Diagram 


Vmax 

58 

ft/s 

Vmin 

19 

ft/s 

Vcruise 

23 

ft/s 

CLmax 

1.3 


n lim 

1.15 


n u it 

1.73 


Pmax (bank angle) 

30 

degrees 

W/S (wing loading) 

.00237 

lb/ft 3 


.0379 

oz/ft 3 

F.S. (factor of safety) 

1.5 
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SECTION VI 


WEIGHTS ANALYSIS 


A. Weights Estimation 


An estimation of component weights was initially made based on standard 
component weight fractions for remote control aircraft of a similar type. 
These provided the initial goal weights for the components which guided the 
design of the aircraft. Based on an overall "target weight" for the aircraft of 3 
pounds (48 ounces), the component weight fractions were determined. These 
are summarized in Table 6-1. 

TABLE 6-1: Initial Estimation of Component Weights 


Component 

Weight (oz) 

Weight Fraction (%) 

Structural Frame 

19.2 

40.0 

Propulsion 

17.3 

36.0 

Avionics 

7.2 

15.0 

Landing Gear 

2.1 

5.0 

Payload 

2.0 

4.0 


After the propulsion system, wing, control surfaces, and fuselage were sized 
and the material selected, a more refined estimation of component weights 
was made. The structural weight (fuselage, empennage, wing) were estimated 
from the structural configuration discussed in Section V. The propulsion 
system analysis resulted in the selection of the motor, batteries and propeller, 
and hence the propulsion weight. Avionics weight was known from the 
knowledge of the servos and speed controller used. The payload was 
designated to be 2 ounces. The landing gear remained a true estimation since 
its design had not been determined. The refined component weights and 
weight fractions are summarized in Table 6-2. 
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TABLE 6-2: Refined Estimation of Component Weights 


Component 

Weight (oz) 

Weight Fraction (%) 

Structures 

21.26 

44.3 

Wing 

1036 

21.6 

Fuselage/ Empennage 

6.4 

13.3 

Skin, Fasteners, Links, etc. 

4.50 

9.4 

Propulsion 

16.19 

33.7 

Propeller 

1.00 

2.0 

Fuel 

5.46 

11.4 

Engine 

9.73 

20.3 

Avionics 

4.15 

8.6 

Landing Gear 

2.40 

5.0 

Payload 

2.00 

4.0 


The avionics components consisted of the receiver, system batteries, and 
servos, while the engine consists of the motor and speed controller. The 
"skin, fasteners, links, etc." component is that material which remained for 
use in construction of the prototype. 


B. Internal Layout 

In determining the internal arrangement of the aircraft, several factors had to 
be considered. Those factors of influence were: 

- effects on center of gravity location 

- proximity requirements 

- simplicity of layout: non-interfering components 

- cooling paths 

- access 

- wing mounting space designation 

With the above considerations, the internal layout was determined. 
However, the locations were able to be changed slightly in the event that 
center of gravity location of the prototype was different from the designed 
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location. That is, there was some flexibility in this layout to accommodate for 
the uncertainties of the construction phase. 

Table 6-3 is a summary of the internal system components, their respective 
functional classifications, and the corresponding sizes and weights: 

TABLE 6-3: Weights and Sizes of Internal Components 



* weight and sizes data is given per battery or servo 
** indicates diameter 


Note that the battery pack is comprised of seven AA batteries, each of which is 
1.95 inches in length and 0.5 inches in diameter. Their combined weight, as 
indicated above, is 5.46 ounces. Also, the system batteries were AAA, 1.73 
inches long and .39 inches in diameter. Note that there are two servos: one 
operates the rudder surface and the other operates the horizontal tail surface. 

To determine the relative location of the internal system components, a 
weight distribution program was utilized (see the computer program listed as 
Item 12 in the Appendix). With fixed values for the structural center of 
gravity, the placement of the internal components was selected. With the 
location of the components already set, the program was used to determine 
the location of the battery pack (whose weight was considered to be ballast). 
As mentioned before, the location of the battery pack and other components 
were allowed to shift slightly so that the overall center of gravity would 
remain in the design location. 
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The design locations of the internal components are summarized in Table 6-4 
and shown in Figure 6-1. 

TABLE 6-4: Design Locations for Internal Components 


Component Center of Gravity Location [in] 

x z y 


motor 

1.13 

2.63 

2.0 

speed controller 

4.0 

2.63 

2.0 

receiver 

6.0 

0.63 

2.0 

system battery 

6.0 

2.28 

2.0 

servo 1 (rudder) 

18.94 

2.44 

1.0 

servo 2 (horiz. tail) 

18.94 

2.44 

3.0 

battery pack 

9.5 

1.5 

2.0 

payload 

9.0 

3.0 

2.0 


x is measured from the nose of the aircraft, increasing toward the tail 
z is measured from the bottom of the fuselage, increasing toward the top 
y is measured from the middle of the fuselage, increasing toward the right 

The servos are placed in the rear so there is no interference with the control 
lines to the rudder and elevator. They will be mounted on 1/16" balsa wood 
strips (servo trays), using wood screws. All other internal components, with 
the exception of the motor and the battery pack, will be similarly mounted. 
To avoid wood splitting at the screw holes, epoxy will be applied in the holes 
prior to inserting the screws. The motor will rest on two triangular, balsa 
wood mounts, and will be secured by elastic straps and screws. The battery 
pack will be secured by Velcro straps, to allow for quick re-location in the 
event that additional modifications are made which affect the center of 
gravity location. Thus, the battery pack provides a means for center of 
gravity travel. 
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Figure 6-1 Internal Layout C onfigurat ion 







C Center of Gravity Estimation 


The placement and location of the center of gravity of each component are 
shown in Figure 6-2. The CG of the overall structure is set at 33% of the chord 
of the wing and the locations of the vertical and horizontal tail leading edges 
are set at 31.7' and 33.2”, respectively. With these values fixed, the placement 
of the internal components was selected as outlined in the previous section. 
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Figure 6-2 
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Center of Gravity Estimation 





SECTION VII 


PERFORMANCE 


A. Lift, Drag and Lift/Drag 


As for every aircraft, the basic size, weight, and cruise velocity of the aircraft 
dictate the lift coefficient necessary to sustain steady, level flight. The 
SCREEM-J4D, with its planform area of 5.46 ft 2 , weight of 3 lbs, and cruise 
velocity of 23 feet per second needs a Cl of 0.874 to maintain steady, level 
flight. The lift coefficient necessary for flight is plotted as a function of 
velocity on Figure 7-1. The drag breakdown described in Section II gives a 
drag polar of Cd = .0319 + Cl 2 / 6827CAR Using the SCREEM-J4D's aspect ratio 
of 11.72 and cruise Cl of .874, the drag coefficient at cruise is then .0564 (see 
Figure 7-2) for a total drag of .194 pounds. Combining this drag and the cruise 
velocity, the power required at cruise is 6.05 Watts. While this is not the 
minimum power required for this aircraft, (which incidentally occurs at 16-18 
feet per second, as can be seen in Figure 7-3) it is only 25% higher than this 
minimum power required. With the lift and drag coefficients thus calculated, 
the maximum lift/drag ratio is 14.0, and as is evident in Figure 7-4, the 
SCREEM-J4D is flying at maximum lift to drag ratio at cruise. 


B. Take-Off and Turning Flight 

The SCREEM-J4D produces a static thrust of 1.186 pounds, which allows for 
takeoff in 45 feet at a takeoff velocity of 23.8 feet per second and C L )max of 1.2. 
The power required to climb at the prescribed rate of 6 feet per second is 37 
Watts, and the relatively high power available of 50 Watts means that the 
SCREEM-J4D will possess the capability to climb at a much higher rate of 15.2 
feet per second. 

In turning flight at the prescribed turn radius of 40 feet, the SCREEM-J4D 
must bank to an angle of 23°, which means that the load factor in the turns 
rises to 1.1. At a roll rate of 28.8 degrees per second, banking to this angle 
requires .8 seconds and 18.4 feet to complete the transition to "fully banked." 
In order to perform this maneuver, the SCREEM-J4D must induce a change 
in angle of attack on the inboard wing of -3° and a change of +3 on the 
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outboard wing. This change is incurred by a maximum rudder deflection of 
25° Reducing the rudder deflection will reduce the roll rate, but even cu ng 
*e "In half will no. severely Inhibit the roll performance of .he 

aircraft. 


C. Endurance and Range 

Given .he battery capacity of 500 mahr the SCREEM-J4D, the maximum 
endu"al« * 221 seconds Broken down into the different sections of the 
i, ; 5 mlfu used in the takeoff phase and 20 mahr is used dmrmg £ 
climb phase, leaving 465 mahr for the steady, level p ase o e 
endurance should be sufficient for the mission, whtch w' jniy requ 
seconds. Combining this endurance with the cru.se velocty o 23 fe P 
second gives a maximum range of 5500 feet, once again well above t 
required range of 1800 feet. Parts B and C of this section are summarized 
Figure 7-5 while Table 7-1 summarizes all of the critical performance data. 

TABLE 7-1 : Critical Performance Data 


Vcruise 
Cl in cruise 
Cd in cruise 
L/Dmax 
Pr in cruise 
Pa 

Pa needed for climb 
R/Cmax 


23 ft/s 
.874 
.0564 

15.5 

6.05 Watts 
50 Watts 
37 Watts 
15 ft/s 


n turn 
p roll rate 

E max 
R max 

Y min glide angle 
R/D 


1.1 
29 “7 sec 

for 6r = 25° 
221 sec 
5500 feet 
4 degrees 

1.6 ft/s 
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Figure 7-5 


Power Profile for Final Design 



Velocity [ ft/s ] 


Motor: 

Propeller: 

Batteries: 


Astro 05 

Tornado 10-6 

7 AA Nickel-Cadmium 


Endurance: 

Range: 

R/C: 

T-0 Distance: 


221 seconds 
5500 feet 
6 ft/s 
45 feet 


D. Landing 


Because the aircraft will be flying at its maximum lift to drag ratio at cruise, its 
glide angle will be maximum at cruise. This maximum glide angle is -4°, 
which will give a minimum rate of descent of 1.6 feet per second. At this rate 
of descent, the descent phase will take approximately 12.5 seconds, requiring 
close to 300 feet from the start of the descent phase until touchdown. This 
distance will be shortened somewhat by descending at a larger angle, with a 
flare upon touchdown to reduce the impact force. Due to the rough surface of 
the runway, it is estimated that for a touchdown on the first half of the 
available runway, the aircraft will come to a complete stop by the end. As a 
precaution, a net will be available after landing if it is evident from flight tests 
that the plane requires too much distance to stop. 
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SECTION VIII 

TECHNOLOGY 

DEMONSTRATOR 



A. Construction Phase 


The product which emerges from the manufacturing phase is rarely identical 
to the product of the design phase. The various problems and solutions 
which were encountered in the construction of the SCREEM-J4D, as well as 
the resulting variations from the design, are presented in this part. 

The aircraft construction was facilitated by the detailed design studies, 
especially regarding the structural framework. The fabrication process was 
aided by actual size blue prints which served as templates for the assembly of 
the fuselage, wing, and empennage. The structure was made of balsa wood, 
spruce and plywood as designed and held together by epoxy and model glue. 
The skin was a MonoKote plastic. 

The weight of the prototype compared very well with the design weight. Each 
component of the prototype weighed approximately as it had been designed 
with the overall weight being .05 pounds less than was designed. This was 
due to weight saving techniques employed during construction. Holes were 
punched in non-critical structural components ~ airfoils, false ribs, trailing 
edge ~ and the empennage was built using fewer braces than designed to save 
weight. The MonoKote plastic was also used on the fuselage to eliminate the 
need for heavier balsa planking. 

The prototype did differ from the design in several ways. First, the center of 
gravity position was designed to be at 33% of the chord. However, the 
prototype had its center of gravity at 40% of the chord. This was due in part to 
the internal arrangement being slightly different. After the first flight test it 
was decided that the CG had to be moved up to the 25% position. This center 
of gravity shift was accomplished by lengthening the nose of the fuselage by 
2.6 inches. This represented a significant difference from the design. 

The chord of the main wing was to be 8.2 inches. However, the blue prints 
were drawn with an 8 inch chord and the trailing edge piece was shortened to 
adjust to the change. 
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The prototype possessed a greater endurance than the design. The design 
batteries were not available for purchase. The batteries that were bought 
differed only in capacity (700 mahr instead of 500 mahr); the voltage and 
weight specifications remained the same. This change only affected the 
endurance, which was significantly increased. 

Another striking difference was related to the wing. The design called for the 
wing to be in three sections with the outboard sections removable for storage. 
The joint which connected the outboards to the main was not considered to 
be structurally sound by an expert model builder late in the construction 
phase. It was suggested that the joint be rebuilt for greater strength. It was at 
this point that the executive decision was made to relax the requirement that 
the aircraft fit in a 2’x2'x4' travel box, which restricted any piece from being 
greater than four feet in length. This enabled AETC to permanently affix the 
outboard sections to the main wing, greatly increasing the wing's strength 
and producing a piece 8 feet in length. 


B. Flight Tests 

The SCREEM-J4D was constantly examined and tested during the 
construction phase. The motor and control surface actuators were checked as 
was the strength of the aircraft's structure. It was decided that before the J4D 
flew its first test flight it should undergo an extensive ground check. 
Described here are the pre-flight test plan and the actual test flight plan. 

Pre-Flight Test Plan : The pre-flight test was intended to check all systems so 

as to insure a successful flight around the required course. Those capabilities 
which were examined were: control surface response; ground roll 
controllability; power production; structural integrity; and lift generation. 

Taxi testing was first performed to see how well the aircraft rolled down the 
runway and was controllable by means of the rear wheel and control surfaces. 
The power production was partially examined by measuring the speed of the 
ground roll. The aircraft was held by the wingtips to test their strength. 
Lastly, the aircraft was allowed to accelerate to take-off velocity and then 
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flown for a few feet very close to the ground to better understand its in-flight 
behavior. 

This test was carefully observed in order to identify problems which should 
be addressed before the actual flight. 

plight Test Plan : The course to be flown is shown in Figures 8-1 and 8-2^ 

Th! course consists of three figure-8 maneuvers preceded by take-off and 
climb and followed by landing within the shaded area. The groun r 
should cover no more than 75 feet with the aircraft becoming airborne by 
point 1. The aircraft then climbs at a 15° angle until it reaches cruise altitude 
of 20 feet at the point marked 2. It then banks left and begins its first figure- 
maneuver at point 3. Point 3a indicates the other end of the i m^neuve 
whose radius of curvature at the ends is 40 feet. After three f.gure-8 s the 
aircraft begins to fly a straight path at point 4. banks around the . end of 
flight zone and begins decent at point 4a, touches down after a flare 
maneuver at point 5, and comes to a complete stop by point 6. 

The zone allowed for this flight is 300 feet long, 120 feet wide, and 25 feet 
high. These dimensions are shown in Figure 8-1 and 8-2. The aci tty use 
encloses a football field and has physical dimensions larger than the allowable 
flight zone. This flight is flown in a controlled atmosphere with limited 

disturbances. 


C. Instructions For Pilot 


The pilot was informed of the capabilities and limitations of the aircraft in 
regards to power, turning capability, bank angles, maximum angle of attack, 
and cruise velocity. He was also told of the flight plan as described above. 


The pilot will position the throttle such that the aircraft will roll at a modest 
speed, but not take off. He will get a feel for the control surfaces as well as for 
the turning capabilities of the rear wheel. After he is satisfied wit e ee , 
he will begin the mission flight. 
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The pilot knows that the maximum take-off distance is 75 feet, which is 
marked on the runway. He will stand anywhere he feels comfortable and 
begin the ground roll at full throttle. If for any reason the aircraft fails to 
achieve flight by the 75-foot mark or it turns out of control, he will abort take- 
off immediately. For the case of a successful take-off, he will climb at full 
throttle to an altitude of 20 feet at which point he will throttle down to 
roughly 50%. When the aircraft is over the end of the runway (even with the 
pylon, at point 2) he will begin a left turn with approximately a 25° bank 
angle. However, if for some reason the aircraft does not reach a 20-foot 
altitude by the pylon, he should begin a turn anyway to avoid leaving the 
allowable flight zone. In this case, he should maintain a full throttle setting 
until he is able to reach 20 feet. If an altitude of 20 feet is unattainable, the 
pilot should land the plane as soon as possible. 

Upon successfully climbing and turning, the pilot will proceed to fly three 
figure-8 maneuvers. The throttle setting for steady, wings-level straight 
portions of flight is approximately 50%. The pilot will increase the power 
slightly in the turns so that the aircraft will not loose altitude in the turns. He 
will fly the straight portions of the flight as accurately as possible so that 
measurements may be taken by the design team. 

If at any time the pilot feels uncomfortable with the qualities of the aircraft — 
i.e. stall likely, dynamically unstable, difficult to maintain altitude, 
unresponsive, ineffective controls, etc. — he is advised to land it as soon as 
possible in any safe location. The condition of the aircraft is of higher priority 
than fulfilling the mission requirements. If the aircraft remains in good 
condition, it may be modified to fly later and successfully meet the 
requirements. 

After completing the three figure-8's, the pilot will fly the aircraft around the 
end, throttle down the motor, descend, and land at the beginning of the 
runway. Upon touchdown, he will completely throttle back the motor and 
guide it to a safe stop. 
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D. Safety Considerations 


The SCREEM-J4D has not passed any official safety test, therefore it -*ould 
not be flown except by an experienced RPV pilot. All spectators shou e a 
safe distance from the flight zone and fully attentive in case of a "runaway 
aircraft. Members of the design team within the flight zone should be aware 
at all times. No one should be at the end of the runway as the aircraft takes 
off or lands or under the aircraft as it makes a turn. The aircraft may sta in 
these circumstances and fall upon someone. 


No one should be within three feet of the aircraft while the motor is » g- 
The propellers are extremely dangerous and often hard to see. o one s ou 
ever touch, pick up, or work on the aircraft unless the battery pack is 

disconnected from the motor. 


During the design flight, NO ONE should be in any position which the pilot 
has designated as unsafe. 

The pilot must realize his restrictions flying indoors; hence, he must take 
protective measures including possibly aborting the mission so as not to 
jeopardize the facility and/or spectators. 

Components inside the aircraft generate heat and may be very hot after 
flight or taxi test. The aircraft should be handled with the idea that 
some components may be of high temperatures. The aircraft should be 

allowed to cool between missions. 
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F. Cost Analysis 

The cost breakdown for the SCREEM-J4D is as follows: 


Item or System 

Propulsion System 

Motor & Gearing 
Propeller 
Batteries 
Speed Controller 

Avionics 

Transmitter 
Receiver 
System Batteries 
Servos 

Actuation Rods 
Hinges and Attachments 
Structures 
Balsa 
Spruce 
Glue 

MonoKote skin 

Labor 

Administrative Costs 

Total _ 


75.00 

2.00 

15.75 

100.00 

200.00 

40.00 

9.00 

88.00 

5.00 

20.00 

40.00 

5.00 

15.00 

20.00 
536.00 

50.00 

$1220.75 


7 AA NiCad 


4 AAAA NiCad 
2 


2 rolls 

160 hours @ $3.75 
copying & binding 
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G. Critical Data Summary 

The final form of the Critical Data Summary is included below. 


DESIGN GOALS: 


EMPENNAGE: 


Re Cruise 

100,000 

Horizontal tail: 


V cruise 

23 ft/s 

Area 

.625 ft 2 

Altitude cruise 

20 ft 

aspect ratio 

2.722 in 

Turn radius 

40 ft 

root chord 

5.75 in 

Endurance 

78s 

tip chord 

5.75 in 

Weight 

48 oz 

taper ratio 

1 



l.e. sweep 

0 

BASIC CONFIG: 


incidence angle 

0 degrees 

Wing Area 

5.46 ft 2 

hor. pos. of 1/4 MAC 

2.9 ft 

Weight 

48 oz 

ver. pos. of 1 /4 MAC 

0 

Wing loading 

8.79 oz/ft 2 

Airfoil section 

flat plate 

length 

37.2 in 



span 

96 in (8 ft) 

Vertical Tail: 


height 

14 in 

Area 

.38 ft 2 



aspect ratio 

1.38 

WING: 

1 (low mount) 

root chord 

7.2 in 

Span 

8 ft 

tip chord 

3.9 in 

Aspect Ratio 

11.72 

taper ratio 

0.54 

Root Chord 

8.2 in 

l.e. sweep 

23.2 degrees 

taper Ratio 

1 

hor. pos. of 1/4 MAC 

2.749 ft 

Airfoil section 

NACA 4415 

Airfoil section 

flat plate 

t/c 

15% 



Incidence angle 

4.5 degrees 

AERODYNAMICS: 


C mac - MAC 

8.2 

Cl max 

1.2 

Hor. pos of 1/4 MAC 

9.017 in 

lift curve slope 

3.72/rad 

Ver. pos of 1/4 MAC 

-4 in 

CDo 

0.0319 

Dihedral 

10 degrees 

efficiency ( e ) 

0.682 

Twist 

0 

Alpha stall 

14.5 (wing) 

Sweep 

0 

Alpha zero lift 

-4 (wing) 



L/Dmax 

15.5 

FUSELAGE: 




Length 

37. in 



Diameter - max 

4 x4 in 



Fineness ratio 

9.3 



Payload volume 

25.17 cu in 
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WEIGHTS: 
Weight total 

48 oz 

C.G. (x-pos) 

.33c 

Avionics 

7.38 

Payload 

2 

Engine 

6.5 

Propeller 

1 

Fuel 

5.46 

Structure 

Wing 

10.36 

Fuselage/emp. 

5.22/1.18 

Landing gear 

2.4 

PROPULSION: 

* 

Type 

Astro 05 

number 

1 

placement 

front 

Pavil max ©engine 

100W / 50 W 

Propeller diameter 

10 in. 

Propeller pitch 

6 in. 

Number of blades 

2 

STABILITY AND 
CONTROL: 
Neutral point 

.514 c 

Static margin %MAC 

18.40% 

Hor tail volume ratio 

0.356 

Vert tail volume ratio 

0.016 

Elevator area 

.208 ft 2 

Elevator max 

15 degrees 

deflection 
Rudder Area 

.21 ft 

Rudder max deflection 

25 degrees 

Aileron Area 

none 

Aileron max deflection 

none 

Cm alpha 

-0.685 [1 /rad] 

Cnbeta 

7.68 E-4 

Cl alpha tail 

3.73 

Cn delta r tail 

5.33 E-4 (neg) 

Cm delta e 

-0.683 [1 /rad] 


PERFORMANCE: 

Vmin 

18.9 ft/s 

Vmax 

58.6 ft/s 

Vstall 

18.9 ft/s 

Range max 

5500 ft 

Endurance Max 

3.9 min 

ROC max 

6.3 ft/s 

Min Glide angle 

4.1 degrees 

T/O distance 

42 ft 

T/O rotation angle 

1 1 degrees 

Landing Distance 

50 ft 

Bank Angle 

30 degrees 

Turn Radius 

28.6 ft 

SYSTEMS: 

Landing gear type 

tail dragger 

Main gear position 

wings 

Main gear length 

5 in 

Main gear tire size 

1 in diameter 

tail gear position 

rudder 

tail gear length 

.5 in 

tail tire size 

.5 in diameter 

engine speed control 

yes 

Control surfaces 

Rudder/elevator 

TECH DEMO: 

Weight total Take off 

49 oz 

Weight empty 

49 oz 

Wing Area 

5.23 ft 2 

Hor. Tail Area 

.625 ft 2 

Vert Tail Area 

.38 ft 2 

C.G. position 

,27 c 

1/4 MAC position 

9.0 in 

Static margin %MAC 

18 % 

V takeoff 

30 ft/s 

Range 

8000 ft 

Endurance 

240 s 

V cruise 

33 ft/s 

Re cruise 

128,000 

Turn radius 

30 ft 
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SECTION IX 


CONCLUSIONS 


A. Disadvantages and Advantages 


While the SCREEM-J4D will accomplish the proposed mission with a margin 
of safety, there are a few short comings to the design. The allowable payload 
is limited to a volume of 25.17 cubic inches and a weight of 2 ounces, to meet 
our design requirements. The J4D has also been designed to meet the specific 
requirements of this proposal. This greatly restricts the design flight envelop 
as compared to general aviation RPVs. In addition to the restriction on the 
flight envelop, the J4D was design to fly indoors, in a controlled 
environment. While this will not restrict the SCREEM-J4D to indoor flights, 
it has not been design to handle hostile environments, including wind gusts 
and adverse weather conditions. This does not preclude the SCREEM-J4D 
from all inclement weather; however, further analysis of its durability is 

suggested. 

The SCREEM-J4D is not without its advantages. Due primarily to the 
constraints imposed on the design by AETC and the initial request for 
proposal, the J4D has been designed as low speed, lightweight, highly 
maneuverable RPV. In addition, it has been design for disassemble and 
transportation within an 2 x 2 x 4 cubic foot volume. A further advantage of 
the disassembly, is the ability to interchange parts. This facilitates the 
replacement of broken or damaged parts. The disassembly will also aid in the 
modification of the J4D for other missions. 


B. Applications 

From the analysis presented in this proposal, the SCREEM-J4D will perform 
the the required station keeping mission within the design constraints as set 
by AETC, with a reliable margin of safety. Due to the stringent constraints, as 
set by AETC, the SCREEM-J4D can be a formidable option for numerous other 
missions. Due to the goals of this mission, the battery pack and thus the 
endurance are held to a minimum. Through an increase in the milliamp- 
hours, corresponding to an increase in battery size, the endurance of the J4D 
can be increased to meet the requirements of a prolonged station keeping 
mission. Unfortunately, the increase in battery size will increase the weight 
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and required cruise velocity. Prolonged station keeping mission may include 
low altitude surveillance or electronic cloaking for military operations. Since 
the SCREEM-J4D has been designed to perform its mission at low speeds, it 
may be flown at speeds comparable to Navy ships. An electronic device may 
then be loaded on board to emit artificial waves that simulate electromagnetic 
emission of a Naval vessel. 

Since the SCREEM-J4D is designed with a non-airbreathing engine, it may be 
easily adapted to work as an atmospheric sampler, in both toxic and non-toxic 
environments. Due to the design cruise velocity, the J4D may be modified for 
function as a high altitude long duration station keeping mission. The 
SCREEM-J4D was initially designed to complete a station keeping mission; 
thus, it has been designed to optimize steady-state circling. The required 
modifications for a prolonged flight mission were discussed above. Recall, 
that these modifications will increase the cruise velocity. The cruise velocity 
will be further increased for applications at high altitudes, due to a decrease in 
density. This decrease in density will correspond to a decrease in lift for a 
constant velocity. Thus, a high altitude mission will require an increase in 
cruise velocity. While it may be necessary to select another engine to meet 
these requirements, the present design only mandates 50% of the design 
engine torque. This short discussion has presented a few of the numerous 
missions that the SCREEM-J4D may be modified to accomplish 


C. Final Word From The Design Team 

This project has progressed rapidly from beginning to end, no doubt to the 
dedication and incessant efforts of the team members. The design has 
evolved from a rough sketch and a lot of hand waving into a polished aircraft 
more than capable of performing the mission for which it was intended. The 
SCREEM-J4D has a singular mission at this time; but its basic design, as we 
have already seen, can be easily modified into one capable of performing a 
number of far-reaching surveillance and reconnaisance missions. Its basic 
airframe is lightweight and durable, lending itself well to effortless 
transportation and on-site assembly. In addition to being extremely quiet, its 
electric motor is perfect for missions in atmospherically hostile 
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environments. We, the members of AETC, support the SCREEM-J4D fully 
and have no doubt that it will be successful in its intended mission and any 
other applications that follow. 



Chris "Duke" Sturgis 
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1 ' > 1 

SALFORD UNIVERSITY FTN77 VER. 228S <USER2>S343728217>STUDENT>DES I GN 1 F77 


"OMP ILER OPTIONS: 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038. 

0039. 

0040. 

0041. 

0042. 

0043. 

0044. 

0045. 

0046. 

0047. 

0048. 

0049. 

0050. 

0051. 

0052. 

0053. 

0054. 

0055. 

0056. 

0057. 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 


LISTING INTL DCLVAR NOMAP CHECK NgBIG LOGL DYNM I NOOFFSET LGO NOAI> 
NOFRN FPN NOLUNFREC NOSILENT N0J3PTIMISE NOIMPURE 


AE 441: AEROSPACE DESIGN 

GROUP D 

PARAMETRIC SWEEP NUMBER 1 
24 JAN 1990 


***** INITIALIZE VALUES ***** 


REAL PIi RE, V, CBAR, RHO, VIS, KINVIS, PHI, G, RADIUS, N 
P I s*4. *ATAN < 1. > 

G»9. 81 


***** SET OUTPUT CODE ***** 


WRITEd,*) 'ENTER OUTPUT CODE: 
READ(1,*> IWR 


FOR SCREEN, "6" FOR PRINTER' 


***** INPUT PRELIMINARY DATA ***** 


WRITEd,*) 'ENSURE THAT ALL QUANTITIES ARE IN METRIC UNITS!! 
WRITEd,*) 'ENTER TARGET REYNOLDS NUMBER 

wr i ?E(iTi> C enter local density' 

WRITE* l^i/^NTER LOCAL FLUID VISCOSITY' 

write! i?i> V 'I nter desired turn radius' 

READ(1,*) RADIUS 


***** PRELIMINARY CALCULATIONS ***** 


KINVIS=VIS/RHO 


..... LDDP THEBAN : C HOR I N I NG ..... 


WRITE (IWR,*) 'REYNOLDS NUMBER = 
DO 20 D=0. , 12. 

RADIUS=6. +D 
WRITE (IWR,*) ' ' 

WRITidWR! *) 'RADIUS = RADIUS* 
WRITE (IWR,*) ' ' 

WRITE( IWR, 85) 

WRITE ( IWR, 87) 

WRITE (IWR,*) ' ' 

DO 10 C=0. , 38. 

CBAR* 8 . 02+( C*. 01 ) 

V s * ( RE*KINVIS ) / CBAR 
PHI=ATAN(V**2. / ( G*RADIUS ) ) 
N-l. /COS (PHI ) 

PHI=PHI*180. /PI 
WRITE( IWR, 86) CBAR, V, PHI, N 
10 CONTINUE 
20 CONTINUE 


***** FORMATS ***** 


85 FORMAT (5X< 'MEAN CHORD', 5X, 'VELOCITY', 5X, 'BANK ANGLE', 5X, 'LOAD FAC 
+0R * ) 

86 FORMAT (7X, F7. 3, 6X, F8. 3, 7X, F7. 3. 0X<F7 3) 

87 FORMAT (8X , 'CMD', 11X, 'CM/SD',8X, CDEG3 ) 


***** STOP AND END ***** 


STOP 

END 


END OF COMPILATION CLOCKED 


. 372 SECONDS 







, _ . _. 1TE.M Z 

c 

Stephen 

T. 

Valentin* (Design Group D) 

( 

1) 

C 


J:._Valentin* 

AE 441: AEROSPACE SYSTEMS DESIGN 
written 02-15-90 

{ 

2) 

C 


( 

3) 

C 


( 

( 

( 

4) 
5 > 
6) 

C 


LIFTING— LINE WING ANALYSIS (Program B) 

C 


LIBRARY 'PLOTLIB' 

( 

7) 




( 

8) 

C 


•**•»*»*«»***»*******» variable declarations ***************** 

( 

9) 



REAL THETA ( 8 ) , THETAP ( 20 > , ALPHA (11). ALPHAa (11). COEFM (8.8) 
REAL SOLN ( 8. 1 1 ) , ALPHMIK8. 1 1 ), CL( 1 1 ) » CDI ( 1 1 ) . Y(20. 1) 

( 

10) 



{ 

11 ) 



< 

12 > 



REAL GAMMA (20. 11 ). LIFT(20. 11 ), DRA0(20. 11 ). CLT(20. 11 ). CDT(20, 11 

{ 

13) 



REAL CD0(20. 11 ), L. LCS, MU 

( 

14) 



INTEGER WK1 ( 8 ) . WK2 (8). NP(ll) 

( 

15) 




< 

16) 

C 


************************ initial input *********************** 

< 

17) 



( 

18) 



WRITE (1,*) 'enter output port number:' 

( 

19) 



READ (1.*) IWR 

( 

20) 



WRITE (1,*) 

( 

21) 

C 


WRITE (1.*) 'enter airfoil zero-lift angle of attack : Cdeg 1 ' 

( 

22) 

C 


READ (l.«) AOL 

( 

23) 




( 

24) 

c 


******************** initial output ************************** 

( 

25) 



( 

26) 



WRITE (IWR,*) ' LIFT DISTRIBUTION DATA' 

( 

27) 



WRITE (IWR,*) 

< 

28) 



WRITE (IWR, 1 > 

( 

29) 

1 


FORMAT (6X, 'ALPHA', 14X, 'CL',12X, 'CDI 12X, 'Lift') 

( 

30) 



WRITE (IWR, 4) 

( 

31 ) 

4 


FORMAT (' Cdeg 1 ', 21 X, ' C 1 b 3 ' > 

( 

32) 



( 

33) 



IOPT— 111 

( 

34) 




< 

35) 

C 


*********************** constants declarations *************** 

< 

36) 




< 

37) 

C 

***** freestream ***** 

( 

38) 



DENSTY-2. 3769E-3 

( 

39) 



VELCTY-23. 0 

( 

40) 



Qinf-O. 5*DENSTY*(VELCTY**2) 

( 

41) 

C 

***** wing ***** 

( 

42) 



CHORD-O. 6833333 

( 

43) 

C 


SPAN-8. 0 

( 

44) 

C 


AR-11. 72 

( 

45) 

C 


WAR E A-SP AN**2 / AR 

( 

46) 



WAREA-5. 5 

( 

47) 



SPAN-WAREA/CHORD 

{ 

48) 

C 

***** airfoil ***** 

( 

49) 

c 


LCS— 5. 02 

( 

50) 



LCS— 5. 5 

< 

51) 



AOL— 4. 0 

( 

52) 




( 

53) 



PI-4. *ATAN( 1. ) 

{ 

54) 



DEGRAD— PI/180. 

( 

55) 



•****••*•»*****»**»**»** initialization of variables ********* 

( 

56) 

c 


( 

57) 




{ 

58) 



MU-CH0RD*LCS/(4. *SPAN) 

( 

59) 



AOL»AOL*DEGRAD 

( 

60) 




( 

61) 



DO 10 1-1,8 

( 

62) 



THETA ( I ) - ( 90. /8. > *REAL ( I ) *DEGRAD 

( 

63) 

10 

CONTINUE 

( 

64) 




( 

65) 



ALPHA ( 1 ) -6. 0*DE0RAD 

( 

66) 



ALPHAa ( 1 ) -ALPHA ( 1 ) -AOL 

( 

67) 



NP ( 1 ) —20 

( 

68) 



DO 20 J-l, 10 

( 

69) 



NP ( J+ 1 ) —20 

( 

70) 



ALPHA ( J+l ) -DEGRAD* ( ALPHA ( J ) /DEORAD+O. 2) 

( 

71) 



ALPHAa ( J+ 1 ) -ALPHA ( J+ 1 ) -AOL 

c 

72) 

20 

CONTINUE 

( 

73) 



********* determine matrix of Fourier coefficients *********** 

( 

74) 

C 


( 

75) 




( 

76) 



DO 40 1-1,8 
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( 

( 

( 

( 

< 

< 

( 

< 

< 

< 

< 

< 

( 

( 

( 

( 

< 

c 

( 

< 

< 

< 

< 

< 

( 

( 

( 

< 

< 

< 

< 

< 

< 

< 

< 

( 

( 

< 

< 

< 

< 

< 

< 

c 

( 

( 

( 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

( 

< 

< 

< 

< 

( 

< 

< 

< 

< 

< 

< 

< 


77) 

78) 

79) 

80) 
81) 
82) 

83) 

84) 
89) 
86 ) 

87) 

88 ) 

89) 

90) 

91) 

92) 

93) 

94) 
99) 

96) 

97) 

98) 

99) 
100 ) 
101 ) 
102 ) 

103) 

104) 
109) 
106) 

107) 

108) 

109) 

110 ) 
111 ) 
112 ) 

113) 

114) 
119) 
116) 

117) 

118) 

119) 

120 ) 
121 ) 
122 ) 

123) 

124) 
129) 
126) 

127) 

128) 

129) 

130) 

131) 

132) 

133) 

134) 
139) 

136) 

137) 

138) 

139) 

140) 

141) 

142) 

143) 

144) 
149) 

146) 

147) 

148) 

149) 

190) 

191 ) 

192) 


30 

40 


90 

60 


70 

80 


DO 30 J-l, 8 

CaeFS(!fb)-<SIN<L*THETA<I)))*(l.+L*riU/SIN(THETA<I))) 

continue 

continue 


DO 60 J-l, 11 

DO 90 1-1,8 . Mtf 

ALPHMU ( I » J ) - ( ALPHAa ( J ) )*MU 
CONTINUE 
continue 


90 

100 

110 


120 


CALL INVRT (COEFM, 8# 8, WK1» WK2) 

CALL MULT (COEFM, ALPHMU, SOLN, 8. 8, 1 1, 8, 8, 8) 


******************* calculate lift parameters ***************** 


DO 80 1-1,11 

C L ( I ) -P I *AR*SOLN ( 1 , I) 

SUM-1. O 
DO 70 J— 2, 8 

SUM-SUM+L* ( SOLN ( J, I ) /SOLN (1,1)) **2 

CONTINUE 

CDKI )»( (CL( I >**2>/(PI*AR) )*SUM 
E-l. /SUM 
CONTINUE 


***************** calculate circulation distribution ********** 


°THETAP ( I ) -REAL ( I ) * ( 90. /20. ) *DE0RAD 
Y(I, 1 ) -COS ( THETAP ( I ) ) 

DO 100 J-l, 11 
SUM-0. 0 
DO 90 K— 1 , 8 

1*2 ♦REAL(K)~1 

OVin? 5 > ^ j > SSCHORD > 

CDoVl, J ) - 1 . 7237 E- 2 - ( 1 . 698 E- 2 ) *CLT ( I, J) + ( 0 . 09021 ) * ( CLT ( I , J 

CDO( I, J)— CDO( I, J>-(0. A2|22)*(CLT(I, J)**3) 

CDO( I, J)-CDO(I, J)+(9. 0229E-2 ) * ( CLT ( I , U)**4> 

CDT( I, J)-CDO( I, J)+(CLT(I, J)**2)/(PI*E*AR) 

DRA0(I« J)*0. 3*DEN8TY*(VELCTY**2)*WAREA*CDT( I, J) 

CONTINUE 

CONTINUE 

CONTINUE 


****************** final output ***************************** 


^WRITE 1 7 iwR?*) ALPHAd ) /DEGRAD. CLd ) , CDI d ), Qinf *WAREA*CL(I ) 
CONTINUE 


130 

140 


WRITE (IWR,*) , , _ . , . 

WRITE (IWR,*) 'Span Efficiency Factor is ,1 
WRITE (IWR,*) . . .. . , , 

WRITE (IWR,*) 'Spanwise Lift Distribution: 

WRITE (IWR,*) 

DO 140 1-1, 11 
WRITE (IWR,*) 

WRITE (IWR,*) 

WRITE ( IWR, * ) 

WRITE (IWR, 2) 

WRITE (IWR,*) 

WRITE (IWR,*) 

Dt ^WRITE' J ( IWR?3) Y( J, 1 ), LIFT( J, I ), DRAQ( J, I ) , CLT ( J, I ) , CDT( J, 1 
CONTINUE 
CONTINUE 


'Incidence-', ALPHAd )/DEGRAD, ' degrees' 
' Clb/ftl C lb 3 


2 

3 


FORMAT ( 2 X, 'Y/(b/2) ', 7X, 'LIFT/unit span',9X, 

• /MM A A # n u //II T / (MV / r f\T f \ 


'DRAG ', 9X, 'CLT', 12X, 'CDT ' ) 


FORMAT (2X, F7. 3, 8X, F7. 3, 8X, F7. 3, 8X, F7. 9, 8X, F7. 3) 
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Valentine (Design Group D) 


2-3 


Stephen 


T. 


< 

( 

( 

( 

( 

< 

< 

( 

C 

< 

( 

( 

( 


153) 

154) C 
153) C 

136) C 

137) 

138) C 

139) 

160) C 

161) C 

162) C 

163) 

164) 
163) 


CALL TP LOT 
CALL TLABEL 
CALL TITLE 

PAUSE 


( I OPT# Y, CLT. 20# 20. 3) 

< 'Y/(b/2> 'Lm Coefficient ' ) 

< 'SPANWISE LOAD DISTRIBUTION (aOL* 


CALL TP LOT ( I OPT# Y, CDT, 20# 20, 3) 

CALL ISABEL < 'Y/(b/2) ', 'Drag Coefficient') 
CALL TITLE ('SPANWISE DRAG DISTRIBUTION (aOL« 


STOP 

END 


4 degrees ) ' ) 


4 degrees ) ' > 



1 1 £>A 3 


•a.*.*##********************************************* 


1 > 
2 ) 

3) 

4) 

5) 

6 ) 

7) 

8 ) 
9) 

10 ) 
1 1 ) 
12 ) 

13) 

14) 

15) 

16) 

17) 

18) 

19) 

20 ) 
21 ) 
22 ) 

23) 

24) 

25) 

26) 

27) 

28) 

29) 

30) 
31 ) 

32) 

33) 

34) 

35) 

36) 

37) 

38) 

39) 

40) 
41 ) 

42) 

43) 

44) 

45) 

46) 

47) 

48) 

49) 

50) 
51 ) 

52) 

53) 

54) 

55) 

56) 

57) 

58) 

59) 

60) 
61 ) 
62) 

63) 

64) 

65) 

66 ) 

67) 

68 ) 

69) 

70) 
71 ) 

72) 

73) 

74) 

75) 

it. \ 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 


************************************************** 
***** LONGITUDINAL STATIC STABILITY ANALYSIS ***** 
***** ERIC T. FICK ***** 

***** 4 MARCH 1990 ***** 

***** DESIGN GROUP D ***** 

************************************************** 


(MING) 
(MING) 
CENTER (MING) 
(MING) 
su stem> 

(MING) 
(MING) 
(MING) 
(MING) 
(MING) 
(MING) 

“MOMENT CURVE SLOPE Cl/rad3 (MING) 


CLOM=LIFT COEFFICIENT AT ZERO ANGLE OF ATTACK 
CLAM=LIFT CURVE SLOPE Cl /rad 3 ^ ACBnnvu . MT . 

CMACM=MOMENT COEFFICIENT ABOUT THE AERODYNAMIC 
XACM=X-POSITION OF THE AERODYNAMIC CENTER Cft3 
{mesured from the origin of the coordinate 
X AC CB= AERODYNAMIC CENTER DIVIDED BY MEAN CHORD 
IH=INC IDENCE ANGLE Cdeg3 
S=PLANFORM AREA Cft**23 
ARM“ASPECT RATO 
CSAR=MEAN CHORD Cft3 
CMOM*ZERO LIFT MOMENT COEFFICIENT 
CMAM 


TO CENTER OF GRAVITY 


IT=INC IDENCE ANGLE Idea 3 
ST “PLANFORM AREA Cft**23 
BT»SPAN C f 1 3 

LT=DISTANCE FROM TAIL A. C. 

ART= ASPECT RATIO 

XACtIx-POSITION D OF THE AERODYNAMIC CENTER Cft3 
{measured from the origin of 'the coor S.i?3 te 
CLOT“LIFT COEFFICIENT AT ZERO ANGLE OF ATTACK 

CLAT=LIFT CURVE SLOPE Cl /rad 3 

CMOT=ZERO LIFT MOMENT COEFFICIENT 
CMAT=MOMENT CURVE SLOPE Cl /rad 3 


(TAIL) 

(TAIL) 

(TAIL) 

(TAIL) 

(TAIL) 

(TAIL) 

(TAIL) 

sgstem> 

(TAIL) 

(TAIL) 

(TAIL) 

(TAIL) 


CMOF“ZERO LIFT MOMENT COEFFICIENT 
CMAF=MOMENT CURVE SLOPE Cl /rad 3 


( FUSELAGE ) 
(FUSELAGE) 


XCGCB “CENTER OF GRAVITY POSITION DIVIDED BY MEAN CHORD 
{measured from the leading edge of the unng> 
RHO=FREESTREAM DENSITY 
VEL=FREESTREAM VELOCITY 
G*FREESTREAM DYNAMIC PRESSURE 
E=DOMNHASH ANGLE _ 

EO=DOWNMASH ANGLE ' AT ZERO ANGLE OF ATTACK A __ A _ iy 
DEO=CHANGE IN DOMNMASH ANGLE MITH ANGLE OF ATTACK 
CMO=MOMENT COEFFICIENT AT ZERO ANGLE OF ATTACK 
CMA=CHANGE IN MOMENT COEFFICIENT MITH ANGLE OF ATTACK 
SFNP“STICK FIXED NEUTRAL POINT/MEAN CHORD 

{measured from the leading edge of the wing> 
SFSM=ST ICK FIXED STATIC MARGIN 

{as a fraction of the mean chord> 

ALFA=TRIM ANGLE Crad3 

ALFD=TRIM ANGLE Cdea3 

CLM“MING LIFT COEFFICIENT 

CLT=TAIL LIFT COEFFICIENT 

CLTOT* TOTAL AIRCRAFT LIFT COEFFICIENT 

LM“MING LIFT GENERATED 

FT“TAIL LIFT GENERATED 

LTOT=TOTAL LIFT GENERATED 

CDO=FORM DRAG COEFFICIENT 

CDTOT=TOTAL DRAG COEFFICIENT 

DTOT=TOTAL DRAG GENERATED 

LODTOT“TOTAL LIFT TO DRAG RATIO 


IN, I T ( 21 ) , CMOM (21, 21 ), VH(21, 21), CM0T<21, 21 ) , ALFA(21, 21 ) 
CMAM (21, 21 ), CMAT(21, 21 ), CMA<21,21 >' ^LFD(21 , 21 ) , CM0(21 , 21 ) 

EO, DE0,ST(21 ), BT(21 ), ART(21 ),CLAT(21), LT(21, 21), CLOT 
CLM(21, 21 ), CLT(21, 21 ), LM<21, 21), FT(21, 21>, LT0T<21, 21 ) 

REAL E(21 , 21 ), CD0(21, 21 ), CLT0T(21, 21 > « CDTOT < 21 , 21 ), DT0T(21, 21 ) 
REAL LODTOT (21 , 21 ), XCGCB (21 ) , SFNP <21 ) , SFSM<21 ) 


REAL 

REAL 

REAL 

REAL 


C- 

C 

C- 


INITIALIZE "HARD" VALUES 


CLOM=. 2596 
CLAW=3 7199 


3-2 


************■*•**.****■*********#*#*.* *********#****.*** 


77) 

78) 

79) 

80) 

81 ) 
82) 

83) 

84) 

85) 

86 ) 

87) 

88 ) 

89) 

90) 

91 ) 

92) 

93) 

94) 

95) 

96) 

97) 

98) 

99) 
100 ) 
101 ) C 
102 ) 

103) 

104) 

105) 

106) 

107) C 

108) 

109) 

110 ) 
111 ) 
112 ) 

1 13) 

114) 

115) 

116) 

117) 

118) 

119) 

120 ) 

121 ) C 

122) C 

123) C 

124) 

125) 

126) 

127) 

128) 

129) 

130) 

131) C 

132) C- 

133) C 

134) 

135) 

136) 

137) C 

138) C- 

139) C 

140) 

141 ) 

142) C 

143) C- 

144) C 

145) 

146) 

147) 

148) 

149) 

150) 

151 ) 
152) 


CMACW— 0. 093 
XACCB-. 25 
S-5. 46 
ARW-1 1. 72 
CBAR-. 68333 
XACT— 2. 935 
CLOT =0. 0 
CBART=. 4792 
CM0F=-0. 0002 
CMAF-O. 0341 
RHO=. 0023769 
P I =4. *ATAN! 1. ) 

1 CONTINUE 

WRITEd,*) 'ENTER OUTPUT CODE: C13 FOR SCREEN, C63 FOR PRINTER' 

READ ( 1 , * ) IWR 

WRITEd,*) 'ENTER WING INCIDENCE ANGLE C DEGREES 3 ' 

READd,*) IW 

WRITEd,*) 'ENTER X-POSITION OF AERODYNAMIC CENTER OF WING CFT3 ' 
WRITEd,*) '{measured from the origin of the coordinate system}' 
READd,*) XACW 

WRITEd,*) 'ENTER FREESTREAM VELOCITY CFT/S3 ' 

READd,*) VEL 
Q=. 5*RH0*VEL**2 


WRITE! IWR, 95) IW 
WRITE( IWR, 96) XACW 
WRITE( IWR, 97) XACT 
WRITE( IWR, 99) VEL 

DO 10 1-1,6 

IT( I ) =-l . 5+REAL ( I - 1 )*. 5 
WRITE! IWR, * ) ' ' 

WRITE! IWR,*) ' ' 

WRITE! IWR, 91 ) IT!I) 

WRITE! IWR, 92) 

WRITE! IWR, *) ' ' 

DO 20 J-l, 21 

XCGC3! J)=. 3333 
ST!J)=. 4+REAL ! J-l ) *. 05 
BT ! J ) -ST ! J ) /CBART 
ART(J)=BT! J)**2. /ST!J) 

CLAT! J)-!2. *P I ) / ! 1 . + !2. /ART! J) ) ) 

'O' VALUES FOR THE MOMENT COEFFICIENT - 

CMOW ! J, I >=CMACW+CLOW* ! XCGCB ! J ) -XACCB > 

LT! J, I > = ! XACT- XACW ) + ! XACCB- XCGCB ( J ) )*CBAR 
VH ! J, I ) — ! ST ! J ) *LT ! J, I ) )/!S*C3AR) 

EO— ! !2. *CLOW > / ! P I * ARW ) >*180. /PI 
DEO— ! 2. *CLAW ) / ! P I * ARW ) 

CMOT ( J, I )— VH ! J, I ) *CLAT ! J ) * (EO+IW— IT! I ) ) *P I /ISO. 
CMO! J, I ) -CMOW ! J, I ) +CMOT ! J, I>+CMOF 

'A' VALUES FOR THE MOMENT COEFFICIENT - 

CMAW ! J, I >=CLAW*( XCGCB! J)-XACCB) 

CMAT ! J, I)— VH!U, I ) *CLAT ! U ) * ! 1 . -DEO) 

CMA ! J, I ) -CMAW ! U, I)+CMAT!U, D+CMAF 


STICK FIXED NEUTRAL POINT 

SFNP ! J ) -XACCB- ! CMAF/CLAW ) +VH! J, I )* ( CLAT! J ) /CLAW > * ! 1 . -DEO) 
SFSM ! J ) -SFNP < J ) -XC GCB ! J ) 


FINDING ANGLE OF ATTACK TO TRIM AND LIFT FORCES- 


ALFA! J, I )=CMO! J, I ) / ! -CMA ! J, I > ) 

ALFD! J, I )-ALFA ! J, I )*180. /PI 

CLWIJ, I ) -CLOW+CLAW* ! ALFA ( J, I ) + ! IW*PI/180. )) 

E! J, I )-(2. *CLW! J, I ) ) / !PI*ARW) 

CLT ! J, I )=CLOT+CLAT!J)*(ALFA! J, I )-( IW*PI/180. )-E!J, 
+ ! IT! I )*PI/180. ) ) 

LW ! J, I ) -CLW ! J, I >*Q*S 
FT! J, I ) -CLT ! J, I )*Q*ST!U) 


I ) 




^^a’*****#*****#*'#*****'*' ■*■**•*■*** *•**#*•# ************** 


153) 

154) 

155) 

156) C 

157) C 

158) C 

159) 

160) 
161 ) 
162) 

163) 

164) C 

165) C 

166) C 

167) 

168) 

169) C 

170) 
171 ) 

172) C 

173) C 

174) C 

175) 

176) 

177) 

178) 

179) 

180) 
181 ) 
182) 

183) 

184) 

185) 

186) 

187) 

188) 

189) 

190) 
191 ) 

192) 

193) 

194) 

195) 

196) 

197) 

198) 

199) 

200 ) 
201 ) 
202 ) 

203) 

204) 

205) 

206) 

207) 

208) 

209) 

210 ) 
211 > 
212 ) 

213) 

214) 

215) 

216) C 

217) C 

218) C 

219) 

220 ) 
221 ) 
222 ) 

223) 

224) 

225) 

226) 

227) 

228) C 


LT0T(J, I)=LM(J, I)+FT(J, I) 

LTOT ( J, I ) =LTOT ( J, I >*. 966 
CLTOT ( U» I)=LTOT(J, I)/(Q*S> 

CALCULATING DRAG FORCES 

CDO( J, I )=(. 1 4885+ ( ST ( J ) *. 008) ) /S 
CDO ( J, I ) =1 . 1 5*CD0( J» I ) 

CDTQT ( U, I ) =CDO ( U. I )+(CLTOT(U, I )**2/(. 6815*PI*ARM) ) 

DTOT ( Ji I ) =CDTOT ( J, I)*Q*S 

LQDTOT ( J, I ) =CLTOT ( J, I >/CDTOT(J, I) 

PRINTING STATEMENTS 

MR I TE ( I MR i 90) ALFD ( J, I >, LTOT(J, I ), VH(U, I ) 

+, LODTOT ( Ji I ), SFNP ( J) , SFSM(U) , CMA< J, I ) > ST ( U ) 

20 CONTINUE 
10 CONTINUE 


PLOTTING STATEMENTS 


2 CONTINUE 

MRITE(1,*> 'ENTER PLOTTING OPTION: C13 ■ LIFT. C23 * DRAG. 

+/D RATIO. C 43 = CM)alfa, C53 = TRIM ANGLE OF ATTACK. C63 * 
♦POINT, C 7 3 = STATIC MARGIN' 

READ ( 1 , * ) NUM 
IF ( NUM. EQ. 1 ) THEN 

CALL TPL0T(-01 1, ST, LTOT, 21. 21, 6) 

CALL TITLE( 'TOTAL LIFT GENERATED Clb3 vs. TAIL AREA Cft**23' 
CALL TLASEL ( 'TAIL AREA Cft**23', 'TOTAL LIFT GENERATED Clb3') 
ELSE IF (NUM. EQ. 2) THEN 
CALL TPLOTC-Ol 1, ST, DTOT, 21, 21, 6) 

CALL TITLE( 'TOTAL DRAG GENERATED Clb3 vs. TAIL AREA Cft**23') 
CALL TLA3EL ( 'TAIL AREA Cft**23', 'TOTAL DRAG GENERATED C 1 b 3 ' ) 
ELSE IF (NUM. EQ. 3) THEN 
CALL TPLOT ( -01 1 , ST, LODTOT, 21, 21, 6) 

CALL TITLE( 'LIFT TO DRAG RATIO vs. TAIL AREA Cft**23'> 

CALL TLA3EL ( 'TAIL AREA Cft**23', 'LIFT TO DRAG RATIO') 

ELSE IF (NUM. EQ. 4) THEN 
CALL TPLOT (-011, ST, CMA 
CALL TITLE ( 'Cm ) a 1 f a vs 
CALL TLABEL ( 'TAIL AREA 
ELSE IF(NUM. EQ. 5) THEN 
CALL TPLOT(-OU, ST, ALFD, 21, 21, 6) „ r4SA 

CALL T I TLE ( 'TR IM ANGLE OF ATTACK Cdeg3 vs. TAIL AREA Cft**23'> 
CALL TLA3EL ( 'TAIL AREA Cft**23', 'TRIM ANGLE OF ATTACK Cd*g3') 
ELSE IF(NUM. EQ. 6) THEN 
CALL TPLOT (-011, ST, SFNP, 21, 21, 1 > 

CALL T I TLE ( 'NEUTRAL POINT vs. TAIL AREA Cft**23'> 

CALL TLABEL ('TAIL AREA Cft**23', 'NEUTRAL POINT') 

ELSE IF(NUM. EQ. 7) THEN 

CALL TPLOT(-OU, ST, SFSM, 21, 21, 1 ) 

CALL TITLE( 'STATIC MARGIN vs. TAIL AREA Cft**23') 

CALL TLABEL( 'TAIL AREA Cft**23'< 'STATIC MARGIN') 


C33 = L 
NEUTRAL 


) 


21 , 21 . 1 ) 

TAIL AREA Cft**23') 
Cft**23 ', 'CM)aIfa ' ) 


ELSE 
END IF 

WRITE(1,*> 'MORE 
READ ( 1 , * > NUMMM 
IF ( NUMMM. EQ. 1) GOTO 2 
MR ITE ( 1 , * ) 'CONTINUE ? 
READ(1,*) NUMM 
IF(NUMM. EQ. 1 ) GOTO 1 


PLOTS ? Cl 3 FOR YES, CNOT 13 FOR NO 


Cl 3 YES, CNOT 13 FOR NO 


90 FORMAT (2X, 
+F8. 4) 

91 FORMAT (2X , 

92 FORMAT (4X, 
+6X, 'SFSM', 

95 FORMAT (2X, 

96 FORMAT (2X, 

97 FORMAT (2X , 
99 FORMAT (2X, 


FORMAT STATEMENTS 

F7. 4, 2X, F6. 4, 2X, F6. 4, 2X, F8. 4, 2X, F8. 4, 2X, F8. 4, 2X, FS. 4, 2X, 


'TAIL INCIDENCE ANGLE IS ', F8. 4. 
'ALFA ' , 4X, 'LIFT ', 3X, 'VolRat',5X, 
5X, 'CM ) a 1 fa ', 4x, 'Stall ' ) 

'MING INCIDENCE ANGLE 
'MING AERODYNAMIC CENTER 
'TAIL AERODYNAMIC CENTER 
'FREESTREAM VELOCITY CFT/SD 


' C DEGREES 3 ' ) 
'LoD ) t ' , 5X , 'SFNP 

', F8. 4) 

', F8. 4) 

', F8. 4) 

', F8. 4) 



********************************* ************* 


1) c 

2) C 
3> C 

4) C 

5) C 

6) 

7) 

8) 

9) C 
10) C 
11 ) 

12 ) 

13) 

14) 

15) 

16) 

17) 

18) C 

19) C 

20) C 
21 ) C 

22) C 

23) C 

24) C 

25) C 

26) C 

27) 

28) 

29) 

30) 

31) C 

32) C 

33) 

34) 

35) 

36) 

37) C 

38) C 

39) C 

40) 

41 ) 

42) 

43) C 

44) C 

45) 

46) C 

47) C 

48) 

49) C 

50) C 
51 ) C 

52) C 

53) 

54) 

55) 20 

56) 10 

57) C 

58) C 

59) 

60) 

61 ) 

62) 

63) 

64) 

65) 

66) 40 

67) 30 

68) C 

69) C 

70) C 
71 ) C 

72) C 

73) 

74) 


***i?*?***£^V******************************* 
SIZ I£° THE VERTICAL TAIL & RUDDER 
™*^ 0R0UP D — JEFFREY UANICIK 


REAL SV(1 10), 
WRITEC1,*) ' 
READ(1,*) IWR 


LV (110), BETA ( 110, 110 
ENTER IWR: ' 


) 


******** fixed parameters ********* 

P 1=4. *ATAN( 1. ) 

SW=5. 46 
B=8. 0 
ARW=11. 72 
ETAV=1 . 0 

CLAV=(2. *PI )/2. 45 
CNBWF=-. 0000721 

********************* ************** 

******* VARYING PARAMETERS ******** 

. ”7 £°^f TRAIN I S ^ S f T BY ®ROUP DESIGN ANO "RULES OF THUMB" 

Iv * i. 5 ft - 2. 0 ft 

Sv » . 38 ft'2 - .60 ft' s 2 

tao = 62 - 72 fiaure 2.20 Nelson 

delta rudder max = 25 degrees 

************************5********** 

WRITEd,*) 'ENTER VALUE OF TAO : ' 

READ(1,*> TAO 

SSSfiM’oSSP VALUE 0F DELTA RUDDER : ' 


DO 10 1=37, 42, 1 

SVC I ) =REAL ( I ) /100. 

DO 20 U=1 8, 20, 1 

LV ( J) =REAL ( J ) / 10. 

***** CALCULATING CNB OF VERTICAL TAIL — eqn. 2. 80 NELSON ***** 

D w;sva?S A vVjt/Iswi B EFFICIEf,CY FAET0R <EMPIRICAL p ORMUL. 

SS^cSvSSoIfSp } J J 8o“ ’ /2 1+0 2+0 009 * ARW 

***** CALCULATING CNB OF ENTIRE AIRPLANE ***** 

CNB=CNBWF+CNBV 

***** CALCULATING CN-DELTA-R — eqn. s 2. 86 & 2. 87 NELSON ***** 
CNDELR=-ETAV*VV*CLAV*TAO*P 1/180. 

***** CALCULATING BETA (THE HEADING ANGLE) ******* 

— I&lf-JS , D0 NE BY SETTING CN » 0 FOR TRIM CONDITION IN 
E 2ySnR N r CN = CNB*BE1 A + CNDELR*DELTAR 
. _ __ BETA ( J, I )=-CNDELR*DELTAR/CNB 
WRITE (IWR,*) CNB, CNDELR 
CONTINUE 
CONTINUE 


WRITE( IWR, *) 'TAO = ' , TAO, 'DELTA RUDDER 
WRITE (IWR,*) ' ' 

WRITE(IWR,*> ' Sv (ft A 2) ', ' Lv 

* ' 


DO 30 1=37, 42, 1 
DO 40 J*18, 20, 1 
WRITEdWR, *) 
CONTINUE 
CONTINUE 


SV( I ), LVCJ), BETA ( J, I) 


************************* 


CALL TPLOT ( -01 1 , 
CALL TITLEC 'BETA 
CALL TLA3ELC 'Sv 


SV, BETA, 100, 100, 6) 
vs. Sv') 

(ft' s 2> ', 'BETA (degrees)' 


', DELTAR 
(ft) ', ' 


STOP 

END 


BETA (DEGREE 



ir£h s’ 


C ********* ***************************** 
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1) c 

2) C 

3) C 

4) C 

5) 

6 ) 

7 ) 

S) C 
9) C 
10 ) 

11) C 

12) C 

13) C 

14) C 

15) 

16) 

17) 

18) 

19) 

20) C 
21 ) 20 
22 ) 10 

23) C 

24) C 

25) 

26) 

27) 

28) C 

29) 

30) 


************************************** 

PROGRAM TO GENERATE Cl vs BETA GRAPH 
— GROUP D JEFFREY L. JANICIK 
***************************************** 

REAL BETA ( 40 ) , CL ( 40, 40 > , GAMMA ( 40 ) 

WRITE(1,*> 'INPUT IWR: ' 

READ(1,*> IWR 

**** ASPECT RATIO EFFECT ON CLB FOUND BY FIGURE 3. 9 **** 

CLBGAM=-. 00027 

**** VARYING PARAMETER — DIHEDRAL ***** 

"RULE OF THUMB" GAMMA * 5 DEGREES - 8 DEGREES 

DO 10 1=1, 10 

GAMMA < I ) =6. +REAL < I ) 

DO 20 J=9, 28, 1 

BET A ( J ) =REAL ( J ) 

CL ( J, I )=CLBGAM*GAMMA( I >*BETA( J) 

WRITE (IWR, *) GAMMA ( I ), 3ETA( J ) , CL < J, I > 

CONTINUE 

CONTINUE 

*«**-»* PLOTTING ROUTINE ******* 

CALL TPLOT ( -01 1 , BETA, CL, 40, 40, 10) 

CALL TITLEC'Cl vs. BETA') 

CALL TLA3EL ( 'BETA (degress)', 'ROLL MOMENT COEFFICIENT Cl') 
************************************** 

STOP 

END 
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EDWARD DIESER 
AERO DESIGN 
TRADE STUDY 
22 MARCH 1990 

REAL IC< 15), JC<25), DIHED( 15), YAWA(25), DAL (25, 15) 

P I =4. *ATAN( 1. ) 

YAWA™ YAW ANGLE 
DAL= DELTA ALPHA 
DIHED=DIHEDRAL ANGLE 

'ENTER IWR (1=TERMINAL, 6*PR INTER ) ' 
READ(1,*) IWR 

DO 10 I»l, 15 
IC ( I ) =6. +REAL( I ) 

DIHED(I)=IC(I )*P 1/180. 


WRITE( IWR, *) 
WRITE (IWR, *) 


'DIHEDRAL ANGLE™ ', IC(I) , 'DEGREES ' 

'YAW ANCLE (DEGREES) DELTA ALPHA (DEGREES)' 


DO 15 J=l , 25 
JC( J)=REAL( J) 

YAWA ( J ) = JC ( J ) *P I / 1 SO. 

DAL ( J, I ) = ATAN ( SIN( YAWA ( J ) )*TAN(DIHED( I ) ) ) 
DAL ( J, I)»DAL(J, I ) * 1 80 . /PI 


15 CONTINUE 

15 WRITE! IWR, *) ' ',JC(J),' ',DAL(J» I) 

10 CONTINUE 


CALL TPL0T(-011, JC, DAL, 25, 25, 15) 

CALL TITLE! 'ANGLE OF ATTACK vs. YAW ANGLE for various WING DIHEDRA 
+L (deg)') 

CALL TLABEL( 'YAW ANGLE (degrees)', 'DELTA ALPHA (degrees)') 

STOP 

END 


EDWARD DIESER 


ITEjM 7 


1) c 

2) C 

3) C 

4) C 
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16) C 

17) C 
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EDWARD DIESER 
AERO DESIGN 
TRADE STUDY 
22 MARCH 1990 

REAL ROLMO< 13), ROLI (13), ALPHA ( 13) , LI (1 3) , LIP (1 3 > , R IP 

G=. 5*. 00237*23. **2 
S=4. / 12. *. 683 
p x =4* AT AN ( 1 . ) 

PC=0. 


WRITE<1,*> 'ENTER IWR ( 1 ^TERMINAL. 6=PR INTER ) ' 
READ (1,*) IWR 

DO 3 K=l, 1 1 
5 ROLMO < K ) =0. 0 
ROLMO ( 1 ) =0. 0 

LI (2)=2 
DO 7 K=3, 13 
7 LI (K)=2+REAL(K-2)*4. 


DO 10 1=1, 13 
ALPHA ( I ) =6. +REAL ( I ) 

CLd >=. 6+ . 1*REAL ( I ) 

WR I TE ( IWR , * ) 'ANGLE OF ATTACK* ', ALPHAC I > , '(deg)' 
WRITEdWR,*) 'PANEL SECTION # LENGTH SECT. 

PC=0. 


DO 20 J=2, 13 
PC=PC+1 

ROLI ( J ) =CL ( I ) *G*S*LI ( J > 

ROLMO ( J ) =ROLMQ ( J- 1 ) +ROL I ( J ) 

LIP < J ) =LI ( J ) /48. 

R IP ( J ) =RQL I ( J ) 

DO 30 K=l, 11 

30 R IP ( J ) =R IP ( J) /ROLMO ( 11) 

20 CONTINUE 

20 WRITEdWR, *) ' ',PC, ' '.LI(J>.' ',ROLI(J) 

WRITEdWR,*) 'ROLL MOMENT*', R0LM0U2) 

10 WRITEdWR, *) ' 

10 CONTINUE 

CALL TPLQT ( —01 1 , ALPHA, ROLMO, 13, 13, 1 ) 

CALL TITLE ( 'ROLL MOMENT vs. ANGLE OF ATTACK') 

CALL TLABELC 'ALPHA (degrees)', 'ROLL MOfCNT (ft. -lb. ) 

STOP 

END 


( 15), CL(15) 


ROLL MOMENT' 


ITEM g 


TORNADO 1 0 — ^ 

Perf ormance Esti mate 

Analysis by si^ie olade eiesi^t theory. 


2 -n-fc) 


DESCRIPTION: 2 Blades; 10 In. Dia.; CUFk •' Airfoil; Blade set at 24.2 Degrees at 2.1 

C L ! 3HT CONDITIONS: Calculations adjusted for neither lacft nor Reynolds r.uibers at 15 *W ■ 
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DATA SHEET 


A) 
B > 
C) 


D > 


E) 


G) 
H > 
I ) 


Propel 1 er Desl gnati on : 
Number of Blades: 2 


TORNADO 10-6 

Diameter: 10 


( Inches ) 


Select one of the following airfoil sections: 

1) INVISCID FLAT PLATE 

2) THIN FLAT PLATE 
3> SYMMETRICAL 

— > 4) CLARK Y 
5) RAF-6 

Blade thickness may be entered as either: 

!> Fraction of chord 
— > 2) Inches 

Blade data may be entered at radial locations speci 
1) Fractional Radius 
— > 2) Inches 


F> Radius* at which blade setting is measured: 


Blade setting (i .e« ref angle for whole blade) : 
Enter the number of radial data positions: (3—9): 


Data Point 
1 


Radi us* 

. 625 

1.125 

n -i c 

«£. • 1 .iJ 

2.625 
3. 125 
3- 625 

4.125 

4.625 


Chord 

. 625 
. 75 
.3125 
. 875 
. 9375 
. 9375 
. 3125 
. 6875 
. 5 


Thi ckness* 

.305 
.225 
. 191 
. 165 
. 147 
. 134 
.114 
.091 
,058 


fied as: 


2. 125 

24.2 


Angle 

56. 8 
40.33 
30. 44 
24.2 

19.99 

16.99 
14.7 
13- 0 
1 1 . 67 


J) Select desired refinment of analysis: 

— > 1) Analysis by simple blade element theory. 

2) Analysis including induced velocity. 

3) Analysis including induced velocity and tip los 

K> These Cl /Cd coefficient adjustments may be selected 
— > 1) No Cl /Cd adjustments 

2) Mach number adjustment 

3) Reynolds number adjustment 

4) Mach and reynolds number adjustments 

L) Select altitude in thousands of feet: .02 


M> Specify one of the following: 

— > 1) Airspeed FIXED at: 

2) Propeller RPM FIXED at: 


15 MPH 


N) Range of Advance Ratio to be used in calculations: 
J min: 0 J max: 1 


-;r*3 1 data point ?u-5t :e les = : L an 3i% radius! at hers nust progress au:^: 
Designation *jst =t«rt *icn a Utter and ^ay "at cental* a cc^via . 
r :r sQua r e bo oUdes I v > bp as last data point. 

-:gU? *i;s: :e spec: usd in degrees, lengths ;n inches. 

♦ 'hits ps as specified U rbes 5 and E. 
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\ T £ M "1 


CASE - GROUPD 
WGT = 3.000 
SREF = 5.460 
RHO = 2.3780E-03 
CLTO = .3000 
CDTO = 4 . 0000E-02 
CLMAX = 1.200 
SMAX = 300.0 
MU = .2000 
DIA = .8333 
BVOLTS = 8.400 
KT = .5500 
KV = 5.0000E-04 
PAPM = 6.0000E-02 
RBAT = 4.7000E-02 
FUSAMP = 20.00 
GEARAT = 2.210 
DT = 5 . 0000E-02 
XMAX = 60.00 
NJ = 10 

0000 .1040 7 . 0000E-02 
2400 9.2000E-02 5.0000E-02 
*2900 8 . 5000E-02 7.9000E-02 
3400 7 . 9000E-02 4.7000E-02 
3900 7 . 2000E-02 4.5000E-02 
4900 5. 4000E-02 3.9000E-02 
5500 4 . 3000E-02 3.5000E-02 
6500 2 . 0000E-02 2.1000E-02 
'7500 -6. 0000E-03 7.0000E-04 
10.00 .0000 - 100.0 
V TAKEOFF = 23.55 
MAX CURRENT DRAW (amps) = 39. 23 
MAX MOTOR POWER (hp) - 
MAX MOTOR POWER (watts) = V 34 - 1 
STATIC THRUST (lb)- 1 - 18 °_ . 7 , 

STATIC CURRENT DRAW (amps)- Id.'J 
STATIC PROP RPS= 99.71 _ 

Stic prop rps- 99.7ip«i- 16 - 731 


V AT TO (FT/SEC) = 23.76 
DISTANCE (FT) =45.06 
BATTERY DRAIN (mahs) = lo.UJ 
ADVANCE RATIO AT TO = .2892 

SSB“iT A TO™EF0^ 3 °R^TRTION) = 1.080 

?Sct“n(LB» AT TO (BEFORE ROTATION) - .3840 
CURRENT DRAW AT TO (AMPS) - 18.03 



GROUP D 
3.0 


! WEIGHT IN LBS 

5.46 


! WING REF AREA IN FT2 

.002378 


! AIR DENSITY SLUG/FT3 

.3 


! CL AT TAKEOFF ATTITUDE 

.04 


! CD AT TAKEOFF ATTITUDE 

1.2 


! CLMAX - 

300. 


! limit on takeoff dist ft 

.2 


! FRICTION COEFFICIENT 

.83333 


! PROP DIAMETER IN FT 

8.4 


! battery voltage 

.55 


! KT IN in-oz/amp 

.0005 


! KV IN volts/rpra 

.06 


! armature resistance 

.047 


! battery resistance 

20. 


! FUSE AMPS - MAX DRAW 

2.21 


! gear ratio 

.05 


! INTEGRATION TIME INCREMENT 

60. 


! limit on take-off time (SEC) 

10 


! number of prop data points 

0 . 

.104 

.07 

.24 

.092 

.05 

.29 

.085 

.079 

.34 

.079 

.047 

.39 

.072 

.045 

.49 

.054 

.039 

.55 

.043 

.035 

.65 

.02 

.021 

.75 

-.006 

.0007 

10. 

0 . 

-100. 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025. 01 

0026. 01 
0027. 01 
0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040. 01 

0041 . 01 

0042. 01 

0043. 01 

0044. 01 
0045 

0046. 

0047. 

0048. 

0049. 

0050 

0051 

0052 

0053 

0054. 

0055. 

0056 

0057 

0058 

0059 


01 

01 

01 

01 

01 


01 

01 

01 

01 


01 

01 


MN T ffl2 u ^C N NOSILlNf C i5 0 ^I?^iI'-N8TSg u g| OFFSET L8 ° N0ANSI N °DED 


rt AEROSPACE STRUCTURES * 

C* Computer program to determine * 

C* the internal load distribution * 

C* on an aircraft wing * 

C* Written by Michael H. Park * 

C ■*•*■*■*■*•♦ -4 

PROORAM STRUC1 

REAL x <30>, CLC50), CD<50), CM ( 50 ) , WC50), PY(50>, PZ(50), DMX(50>, 
+ VYC50), VZ(50), MX (50), MY (50), MZ(50). C ( 50 ) 


10 


20 


30 


90 

1 


S£?nroa IT 7™' EILEs'WHSJGI. DATA'. STATUS* 'OLD ' ) 
R|AD(99, *)R0, ALPHA. V. CR. CT. N 

CLOSE TunItI99 ) ’ ' CL ‘ I ’ ‘ CO < I > ■ CM < I > . M 1 1 ) . 1 = 1 , N ) 


Q=0. 5-»R0*V*-»2 
P I =4. 0*ATAN( 1. 0) 
A=ALPHA*P 1/130. 0 


WRITE! 1, *) 'ENTER 
READ ( 1, * ) IWR 


IWR: 


DO 


5 1 = 1, N 

C ( I >=CR-( (CR-C 
CONTINUE 


T)/X(N))*X(I) 


D0 


10. 1 = 1, N 

p X\\ C ( I >*!CL! I >*COS(A>+CD( I )*SIN!A) >-W! I )*COS(A) 
nMY? tT 2^ C A Jit CCL( 1 >*SIN<A)-CD( I )*COS(A) ) -W ! I ) *S IN ( A ) 

CONTINUE _Q * C ( 1 °* CM( I )-CL! I )*C0S!A>/4. O-CD ( I ) -»SIN ( A ) / 4 


0) 


VY ( N ) =0. 
VZ CN)=0 
MX ( N ) *0. 
MY ! N ) =0. 
MZ(N)*0. 


DO 20 r=N-l, 1, -1 

ri*w-r5 J + H + J£X ( J + 1 > + PY!I))*!X! I+1J-X! I) )/2. 
MYf T + J -*-(PZC 1 + 1 )+PZ(I ) )*<X( 1 + 1 )-X< I ) )/2. 

CONTINUE *** 11 +(DMX ( 1 + 1 )+DMX <I))*<X(I + 1)-X(I>)/ 


DO 30, I =N— 1 , 1, -1 

Mwf! = tl^J + 1) + (VY{I + 1)+VY(I) >* ( X(I + l)-X(I))/2 0 
CONTINUE =MY( 1+1 +<VZ( 1+1 >+VZ( 1 ) )*<X( 1+1 >-X< I ) )/2. 0 


I! Li L.jti itLiiziiii 


Mr ' 


DO 90 1=1, N 

CONTI Nu| ( 1 WR ' 1 > x < 1 > ' VY ( I ) , VZ ( I ) , MX C I ) , MY ( I ) , MZ ( I ) 
FORMAT! IX, F5. 2, IX, F9. 5, IX, F9. 5, IX, F9. 5, IX, F9. 5, IX, F9. 5) 


STOP 

END 


OF COMPILATION CLOCKED 


• 572 SECONDS 



X 

V<j ( lbs ) 

Vz 

0. 00 

16. 23663 

1. 29189 

0. 31 

14. 83717 

1. 18189 

0. 63 

13. 44794 

1. 07263 

0. 93 

12. 07876 

0. 96484 

1. 24 

10. 74015 

0. 85929 

1. S3 

9. 44182 

0. 73669 

1 82 

8. 19319 

0. 63776 

2 09 

7. 00408 

0. 36321 

2 35 

3. 88393 

0. 47376 

2 60 

4. 84188 

0. 39013 

2 83 

3. 88608 

0. 31299 

3 04 

3. 02453 

0. 24303 

3. 24 

2. 26480 

0 18096 

3. 41 

1. 61361 

0. 12748 

3. 56 

1. 07557 

0. 08321 

3 70 

0. 65281 

0. 04858 

3. 30 

0. 34374 

0 . 0^3 / 2 

3. 89 

0. 14259 

0 . 00831 

3. 95 

0. 03600 

0 00119 

3. 99 

0. 00000 

o ocooo 


Mx ( f t-lbs > 

My 

Mz 

-20. 69218 

2. 41726 

30. 

-18. 81430 

2. 02908 

25. 

-16. 95127 

1. 67731 

20. 

-15. 11737 

1. 36366 

17. 

-13. 32781 

1 . 08796 

13. 

-11. 59669 

0. 84989 

10. 

-9. 93779 

0. 64816 

8. 

-8. 36353 

0. 48087 

6. 

-6. 89396 

0. 34346 

4. 

-5. 33662 

0. 23893 

2. 

-4. 30560 

0. 15784 

1 . 

-3. 21270 

0. 09857 

1 . 

-2. 26930 

0. 05735 

0. 

-1. 48536 

0. 03044 

0. 

-0. 86707 

0. 01427 

0. 

-0. 41562 

0. 00561 

0. 

-0. 12537 

0. 00166 

0. 

0. 01781 

0. 00031 

0. 

0. 04280 

0. 00002 

0. 

0. 00000 

0. 00000 

0. 


24091 

36481 

93403 

02190 

37303 

39984 

08474 

00241 

31949 

99679 

99028 

25361 

73938 

40103 

19466 

08104 

02686 

00614 

00066 

00000 



0 . 00237 , 7 . 6 , 23 . 0 , 8 . 


«=. / 


3 . 2 , 20 


/ ^ - J> 


( 

C 

( 

C 

( 

( 

( 

( 

< 

c 

( 

( 

{ 

{ 

( 

( 

( 

( 

c 

< 

( 


1 ) 0. 00237, 7. 6, 23 

2) 0. 00000, 0. 94210 

3) 0.31384,0.94164 

4) 0.62572,0.94020, 

5) 0.93380,0.93769, 

6) 1.23608,0.93406, 

7) 1.53072,0.92920, 

8) 1.81596,0.92279, 

9) 2.09000,0.91434, 

10) 2. 351 16, 0. 90336, 

11) 2.59780,0.88944, 

12) 2. 82844, 0. 87191, 

13) 3.04164,0.84938, 

14) 3.23608,0.81984, 

15) 3. 41056, 0. 78121, 

16) 3. 56404, 0. 73152, 

17) 3. 69552, 0 66802, 
13) 3.80424,0.58601, 

19) 3.88948,0.47931, 

20) 3.95076,0.34342, 
21 ) 3. 98768, 0 . 13015, 


), 8 . 2 . 3 . 2 , 20 
0. 04717, -0. 093, 0. 375 
0. 04712, -0. 093, 0. 375 
0. 04699, -0. 093, 0. 375 
0. 04677, -0. 093, 0. 375 
0. 04644, -0. 093, 0. 375 
0. 04602, -0. 093, 0. 375 
0. 04547, -0. 093, 0. 375 
0. 04476, -0. 093, 0. 375 
0. 04386, -0. 093, 0. 375 
0. 04278, -0. 093, 0. 375 
0. 04148, -0. 093, 0. 375 
0. 03992, -0. 093, 0. 375 
0. 03803, -0. 093, 0. 375 
0. 03582, -0. 093, 0. 375 
0. 03332, -0. 093, 0. 375 
0. 03057, -0. 093, 0. 375 
0. 02754, -0. 093, 0. 375 
0 02417, -0. 093, 0. 375 
0. 02049, -0. C93, 0. 375 
0. 01727, -0. 093, 0. 375 



i-FORD UNIVERSITY FTN77 VER 

MPILER OPTIONS. 


I T E M i i 

228b -:USER2;>S3 1 0768966>STUDENT>STR0Cij. K77 


20 : 1 


DCLVAR NOMAP CHECH NOBIG LOGL DYNM NOOFFSET LGO NOAh 
NOFRN FPN NOLUNFREC NOS I LENT NO_OPTIMISE NO IMPURE 


0001 

0002 

0003 

0004 
0003 
0006 

0007 

0008 

0009 

0010 
001 1 
0012 
0013 
C014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034. 01 

0035. 01 

0036. 01 
0037 
0038. 01 

0039. 01 

0040. 01 

0041. 01 

0042. 01 

0043 

0044 

0045 

0046 

0047. 01 

0048. 01 

0049. 01 

0050. 01 

0051. 01 

0052. 01 

0053. 01 

0054. 01 

0055. 01 
0056 

0057. 01 

0058. 01 

0059. 01 

0060. 01 

0061 . 01 
0062. 01 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 
C074 
0075 


C 444 4 4444 444*. ■* 4 * * * ^ * * *** * * * *. + # # * j,. # * 

C* AE446 - AEROSPACE STRUCTURES * 

C* Co/np utsr program to determine * 

C* the centroid location and sec- 
C* tion properties relative to the 
C* centroid for a nouniform airfoil 
C* and the direct stress due to 
C* axial* bending and thermal loads 
C-» Written by Michael H. Park 

C *■*•*■»*■» ■»■»■»•*•»■»■»•»****»»**■*•»■»•» •**■*•****•*•*♦** 

PROGRAM STRUC4 

rk,ii. 50) ' A (50). YC50), Z (50), PT(50), MZT(50), 
+ , SIG(50), ALPHC50), MY, MZ 

FILE> 'PRDPS5. DATA', STATUS* 'OLD ' ) 

READ (99, * )N 

F£ A£i.? 9 ' * > < A( I ) , Z C I , Y ( I ) , E ( I ) , ALPH ( I ) , I = 1 , N 
CLOSE (UNIT=99 ) 


* 

* 

* 

* 

* 

* 

4 

* 


MYT ( 50 ) , IYY, IZZ, IYZ 
) 


1 


WRITE(1, 4) 'ENTER IWR: ' 

READ ( 1,^) IWR 
WRITE( 1, *) 'ENTER P, My, 
READ ( 1, 4>p, MY, MZ, DT 


Mz, 


DT: ' 


ATCT =0. 0 
YA=0. 0 
ZA=0. 0 
I YY=0. 0 
IZZ=0. o 
IYZ=0. 0 


1 = 1 

EREF=E( I ) 

DO 5 1=1, N-l 

„ IF <E( I ). LT. E( 1 + 1 > )EREF=E( I ) 

5 CONTINUE 


DO 10, 1 = 1, N 

ATQT=ATOT+ (Ed) /EREF ) -*A ( I ) 
VA=YA+ ( E ( I ) /EREF ) *Y ( I > *A ( I ) 
ZA=ZA+(E( I ) /EREF ) *Z ( I ) *A ( I ) 
CON i INUE 

Y3AR=YA/AT0T 

Z3AR=ZA/AT0T 


20 


DO 20, 1 = 1, N 

IYY=IYY+(E( I ) /EREF ) *A ( I ) * ( Z ( I ) -ZBAR ) 442 
I ZZ=IZZ+(E( I ) /EREF )*A(I)*(Y(I ) -YBAR ) 442 
I/Z=IYZ+(E(I) /EREF ) 4 A ( I )*( Z ( I ) -ZBAR ) ■» ( Y ( I ) -YBAR ) 

PTC I )=E( I ) 4ALPH ( I >*DT 
NYT ( I ) =E ( I ) *A ( I ) *Z ( I ) *DT 
MZT ( I ) =E ( I ) *A ( I ) 4Y ( I ) *DT 
CONTINUE 


30 


DO 30, 1=1 , N 

B=E( I ) * (P+PT ( I ) ) /(EREF*ATOT) 

C=(E( I ) /EREF ) 4 ( ( MZ— MZT ( I ) ) 4 I YY+ ( MY+MYT ( I) )*IYZ)*Y< 
D=(E( I )/EREF)4( ( MY+MYT ( I ) ) * I ZZ + ( MZ-MZT ( I ) >*IYZ)*Z( 
SIG( I ) =3+ ( D-C ) / ( IYY4I ZZ-IYZ**2)-E( I )4ALPH( I )4DT 


I ) 
I ) 


WRITEC IWR, 4) 'YBAR 
WRI TE ( IWR, 4) 'ZBAR 
WR I TE ( I WR , 4 ) ' I y y 
WRITE (IWR, 4) 'l£l 
WR I TE ( I WR , 4 ) ' I u i 
WR I TE (IWR, 4) 'P= 
WR I TE ( IWR, 4) 'My = 
WRITEC IWR, 4) 'Mz = 
WRITE(IWR, 4) 'DT = 


= '- Y3AR, 'in. 

= ' , Z3 AR, 'in. 
', IYY, ' in~4 ' 
', IZZ, ' i n~4 ' 
', IYZ, ' i n' v 4 ' 
, P, 'lbs ' 

', MY, 'in-lbs ' 
', MZ, 'in-lbs ' 
' , DT , 'degrees 


F ' 


WRITEdWR, 4) ' ' 

WPITEC IWR, 98) 

WR I TE ( IWR, 4) ' — 


SALFORD UNIVERSITY FTN77 VER. 228S 


ITEM \X 

<USER2>S343728217>STUDENT?CG2. P77 20: O: 


COMPILER OPTIONS: 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056. 01 

0057. 01 

0058. 01 

0059. 01 
0060 
0061 
0062 
0063 

0064. 01 

0065. 01 
0066 

0067 

0068 
0069 

0070. 01 

0071. 01 

0072. 01 

0073 

0074 

0075 


C 

C 

C 

C 

CC 

c 

c 


20 


bniaM N it>J. N Ik DCi -VAR NOMAP CHECK NOBIO LOGL DYNM NOOFFSET LGO NC 
NOFRN FPN NOLUNFREC NOSILENT NO OPTIMISE NOIMPURE 


REAL IY, X!30), Z!30) , WOO) 


READ?i 1 iT > IWR NTER 0UTPUT C0DE: C13 for SCREEN, C63 FOR PRINTER 

continue 

XC=0 

ZC=0. 

WT=0. 

I Y=0. 


125 

625 

5 

0 

625 

23 

0 

625 

95 

0 

28 

0 


0 

0 

0 


X(1 > = 1. 

Z! 1 )=2. 

W< 1 ) =7. 

X ( 2 ) =4. 

Z(2 )=2. 

W ( 2 ) =3. 

X!3 >=6. 

Z (3 )=0. 

W ! 3 ) =0. 

X ( 4 ) =6. 

Z ( 4 ) =2. 

W ( 4 ) =2. 

X ( 5 ) = 18. 94 
Z ( 5 > =2. 44 
W(5)=l. 

X!6 )=9. 

Z ! 6 ) =3. 

W!6 >=2. 

X ( 7 ) =14. 834 
Z ( 7 ) =2. 25 
W C 7 ) =5. 22 
X (8 > = 12. 62 
Z(8)=0. 75 
W<8> = 10. 36 
X ( 9 ) =33. 67 
Z ( 9 ) =4. 25 
W ( 9 ) =0. 4574 
X( 10) =35. 58 
Z( 10)=4. 0 
W( 10)=0. 7226 
X ! 1 1 )=9. 0 
Z( 1 1 )=-3. 0 
W ( 1 1 ) =2. 4 
W< 1 2 ) =5. 46 
WRITE! 1, *) 'ENTER 
READ! 1, *)X( 12), Z 
N= 1 2 
DO 20 J=l, N 

XC=XC+X!U)*W! J) 
ZC=ZC+Z! J)*W! J) 
WT=WT+W! J) 
CONTINUE 
XCG=XC/WT 
ZCG=ZC/WT 


X- AND 
! 12 ) 


Z-POSITIONS OF BATTERY PACK 


DO 30 K=1 , N 

I Y=I Y+ ! ! X ( K ) -XCG ) **2+ ! Z ! K ) -ZCG ) **2 )*W!K) 
30 CONTINUE 

IY=IY/!144. *16. *32. 2) 

WRITE! IWR, 90) 


DO 40 L=l, N 

„ WRITE! IWR, 91 ) L, X ( L ) , Z ! L ) , W ! L ) 
40 CONTINUE 


WR I TE ! I WR , * ) ' ' 

WRITE! IWR, 92) XCG 
WR ITE! IWR, 93) ZCG 


I 2 .- 2 - 


0076 


WR I TE ( I WR » 95 ) WT 

0077 


WRITE( IWR, 94) IY 

0070 


WRlTEtl,#) 'CONTINUE ? C 1 3 FOR YES' 

0079 


READ ( 1, *) NN 

0080 


IF ( NN. EQ. 1 ) GOTO 1 

0081 

90 

FORMAT (2X* 'COMPONENT '• 3X, 'X-P0S',4X, ' 

0082 

91 

FORMAT < 5X# 12, 6X, F7. 4, 2X, F7. 4, 2X, F8. 4) 

0083 

92 

FORMAT (2X, 'X-POSITION OF THE C. G. = 

0084 

93 

FORMAT (2X, 'Z-POSITION of THE c.g. 

0085 

94 

FORMAT ( 2X , 'MOMENT A30UT Y-AXIS <Iy> =* 

0086 

0087 

95 

FORMAT ( 2X, 'TOTAL WEIGHT * 

0088 


STOP 

0089 


END 


2X« 


F7. 
F7. 
F7. 4, 
F7. 4, 


4, 

4, 


IX 

IX 

IX 

IX 


'WEIGHT CozI ' ) 

, 'INCHES' ) 

, 'INCHES' ) 

, 'SLUG*FT**2 ' ) 
, 'OZ ' > 


END OF COMPILATION CLOCKED 


. 500 SECONDS 



COMPONENT 

X-PQS 

Z-PCS 

WEIGHT Co 

1 

1. 1250 

2. 6250 

7. 5000 

2 

4. 0000 

2. 6250 

3. 2300 

3 

6. 0000 

0. 6250 

0. 9500 

4 

6. 0000 

2. 2800 

2. 0000 

5 

18. 9400 

2. 4400 

1. 2000 

6 

9. 0000 

3. 0000 

2. 0000 

7 

14. 8340 

2. 2500 

5. 2200 

8 

12. 6200 

0. 7500 

10. 3600 

9 

33. 6700 

4. 2500 

0. 4574 

10 

35. 5800 

4. 0000 

0. 7226 

11 

9. 0000 

-3. 0000 

2. 4000 

12 

9. 5000 

1. 5000 

5. 4600 


X-POSITION OF THE C. 0. 
Z-POSITION OF THE C. G. 
TOTAL WEIGHT 
MOMENT A30UT Y-AXIS <Iy> 


9. 6998 INCHES 
1.6286 INCHES 
41. 5000 OZ 
0. 0251 SLUG*FT**2 



NOMENCLATURE 


AETC Aerospace Engineers for the Twenty-First Century 
AR Aspect Ratio 

a Angle of Attack 

b Wing Span 

(3 Propeller Blade Angle 

Yaw Angle 

c Wing Mean Chord 

Cd Drag Coefficient 

Cdo Parasite Drag Coefficient 

CG Center of Gravity 

Cl Lift Coefficient 

Ci Roll Moment Coefficient 

CL a t Change in Lift Coefficient with Angle of Attack for the Tail 

CL av Change in Lift Coefficient with Angle of Attack for the Vertical Tail 

CLq Lift Coefficient at Zero Angle of Attack 
CLa Change in Lift Coefficient with Angle of Attack 
Cjp Roll Moment due to Yaw 

CL W Wing Lift Coefficient 

C m Pitch Moment Coefficient 

Cm acw Moment Coefficient About the Aerodynamic Center of the Wing 
Cm a Change in Moment Coefficient with Angle of Attack 

Cm a f Change in Moment Coefficient for the Fuselage 

Cm a t Change in Moment Coefficient for the Tail 

Cm tt w Change in Moment Coefficient for the Wing 

Cm Q Total Zero Lift Moment Coefficient 

Cm c f Zero Lift Moment Coefficient for the Fuselage 

Cm G t Zero Lift Moment Coefficient for the Tail 

Cmow Zero Lift Moment Coefficient for the Wing 

Cmge Change in Moment Coefficient with Rudder Deflection 

Cn Yaw Moment Coefficient 

Cn|3 Yaw Moment Coefficient due to Yaw Angle b 
Cng r Rudder Control Power 

Cp Propeller Power Coefficient 



Q Propeller Thrust Coefficient 

e efficiency factor 

£o Downwash Angle at Zero Angle of Attack 

P Dihedral 

it Tail Incidence Angle 

iw Wing Incidence Angle 

J Advance Ratio 

J4D SCREEM-J4D aircraft 

lh Horizontal Tail Moment Arm 

l v Vertical Tail Moment Arm 

mahr Milliamp-hours 

n Load Factor 

ti Efficiency Factor 

PA Power Available 

PR Power Required 

Re Reynolds Number 

SFSM Stick Fixed Static Margin 

Sh Horizontal Tail Surface Area 

S v Vertical Tail Surface Area 

Sw Wing Surface Area 

X Flap Effectiveness 

Vh Volume Ratio of Horizontal T ail 

Vy Volume Ratio of Vertical Tail 

W Weight 

X ac Distance From Leading Edge to Aerodynamic Center of Wing 
Xqg Distance From Wing Aerodynamic Center to Center of Gravity 
Xnp Stick Fixed Neutral Point 



